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ABSTRACT 

We use the conditional luminosity function $(L|M)dL, which gives the number of 
galaxies with luminosities in the range L±dL/2 that reside in a halo of mass M, to link 
the distribution of galaxies to that of dark matter haloes. Starting from the number 
density of dark matter haloes predicted by current models of structure formation, 
we seek the form of "I>(_L|M) that reproduces the galaxy luminosity function and the 
luminosity dependence of the galaxy clustering strength. We test the models of $(L|Af) 
by comparing the resulting mass-to-light ratios with constraints from the TuUy-Fisher 
(TF) relation and from galaxy clusters. A comparison between model predictions and 
current observations yields a number of stringent constraints on both galaxy formation 
and cosmology. In particular, this method can break the degeneracy between flo and 
the power-spectrum normalization erg, inherent in current weak-lensing and cluster- 
abundance studies. For flat ACDM cosmogonies with erg normalized by recent weak 
gravitational lensing observations, the best results are obtained for flo ^ 0.3; flo ^ 0.2 
leads to too large galaxy correlation lengths, while Hq 0.4 gives too high mass- 
to-light ratios to match the observed TF relation. The best-fit model for the ACDM 
concordance cosmology with — 0.3 and 17a = 0.7 predicts mass-to-light ratios that 
are slightly too high to match the TF relation. We discuss a number of possible effects 
that might remedy this problem, such as small modifications of erg and the Hubble 
parameter with respect to the concordance values, the assumption that the universe 
is dominated by warm dark matter, systematic errors in current observational data, 
and the existence of dark galaxies. We use the conditional luminosity function derived 
from the present data to predict several statistics about the distribution of galaxy 
light in the local Universe. We show that roughly 50 percent of all light is produced in 
haloes less massive than 2 x 10^^ H^^Mq. We also derive the probability distribution 
P{M\L)dM that a galaxy of luminosity L resides in a halo with virial masses in the 
range M ± dM/2. 

Key words: galaxies: formation — galaxies: clusters — large-scale structures: cos- 
mology: theory — dark matter 



1 INTRODUCTION 

It is well established that galaxies reside in extended dark 
matter haloes. Because of this, the standard assumptions 
in the current paradigm of structure formation are that a 
weakly interacting mass component (such as the cold dark 
matter, hereafter CDM) dominates the mass of the Universe 
and that galaxies form by the cooling and condensation of 
gas in dark matter haloes. Numerical simulations and ana- 
lytical models can now give us a fairly detailed picture of 
the abundance and clustering of CDM haloes (e.g. Mo & 
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White 2002). However, in order to build a coherent picture 
of galaxy formation, we need to be able to link these proper- 
ties of the halo population to those of the galaxy population. 
In other words, we need to know how dark matter haloes are 
populated by galaxies. 

Three different approaches have been used in the past 
to link galaxies to dark matter haloes. The first method tries 
to infer the properties of a halo from the properties of its 
galaxy population. The second method relies on ab initio 
models for the formation of galaxies, and the third method 
attempts to link the distributions of dark matter haloes and 
galaxies using a more statistical approach. 

Mass estimates of individual dark matter haloes can be 
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obtained from galaxy kinematics, from gravitational lens- 
ing, or from studying the X-ray haloes. However, each of 
these methods has its own shortcomings: kinematics gen- 
erally only probe the inner regions of dark matter haloes, 
results from gravitational lensing depend rather sensitively 
on the models adopted for the lensing systems, and X-ray 
measurements have to rely on the assumption that the gas 
is in hydrostatic equilibrium. Although significant progress 
has been made in each of these areas, the main limitation 
is that all these methods yield mass-to-light ratios for in- 
dividual objects only. Because of this, various authors have 
attempted a more global approach by using a combination 
of galaxy luminosity functions and luminosity-velocity re- 
lations (such as the TuUy-Fisher relation) to construct a 
distribution function for halo circular velocities (or masses) , 
which in principle can be compared directly to the results 
from A'^-body simulations (i.e., Shimasaku 1993; Kauffmami 
et al. 1993; Newman & Davis 2000; Gonzales et al. 2000; 
Bullock et al. 2001a; Kochanek 2001). Such analyses can be 
used to put constraints on the typical mass-to-light ratio for 
a population of galaxies, such as spiral galaxies. However, 
this method needs to make a rather unreliable conversion 
from an observed velocity measure, such as the maximum 
rotation velocity of a disk galaxy, to the virial velocity of the 
dark matter halo (see discussions in Gonzales et al. 2000; 
Kochanek 2001; van den Bosch 2002). 

An alternative approach is to build ab initio models for 
the formation of galaxies, using either large numerical sim- 
ulations (e.g., Katz, Weinberg & Hernquist 1996; Fardal et 
al. 2001; Pearce et al. 2000; Kay et al. 2001), (semi-) an- 
alytical models (e.g.. White & Rees 1978; White & Frenk 
1991; Kauffmann, White & Guiderdoni 1993; Somerville & 
Primack 1999; Cole et al. 2000; Benson et al. 2002; van den 
Bosch 2001, 2002), or a combination of both (Kauffmann, 
Nusser & Steinmetz 1997; Kauffmann et al. 1999; Diaferio 
et al. 1999; Benson et al. 2000, 2001, Springel et al. 2001; 
Mathis et al. 2001). Although these methods have yielded 
many useful predictions about how dark matter haloes are 
populated by galaxies of different intrinsic properties, many 
of the physical processes (such as star formation and feed- 
back) involved are still poorly understood at the present 
and so some uncertain assumptions have to be made about 
a number of model ingredients such as the efficiency of star 
formation and feedback, the stellar initial mass function, and 
the amount of dust extinction. 

With the advent of large, homogeneous data sets such as 
the Two Micron All Sky Survey (2MASS), the Two Degree 
Field Galaxy Redshift Survey (2dFGRS), and the Sloan Dig- 
ital Sky Survey (SDSS), it becomes possible to address the 
link between dark matter haloes and galaxies using a more 
statistical approach. The population of galaxies in individ- 
ual haloes can formally be represented by a halo-occupation 
function, P{N\A{[), which gives the probability that a halo 
of mass M contains N galaxies (with certain intrinsic prop- 
erties). Together with an assumption about the spatial dis- 
tribution of galaxies in individual haloes, P{N\M) is suf- 
ficient to transform the abundance and clustering of the 
halo population into those of the galaxy population. Nu- 
merous studies in the past have used models of P{N\M) to 
study the bias and various other clustering statistics of the 
galaxy distribution (e.g., Jing, Mo & Borner 1998; Peacock 
& Smith 2000; Seljak 2000; Scoccimarro et al. 2001; White 



2001; Jing, Borner & Suto 2002; Berlind & Weinberg 2002; 
Bullock, Wechsler & Somerville 2002; Scranton 2002; Kang 
et al. 2002; Marinoni & Hudson 2002). The disadvantage of 
this method is that it does not provide a physical under- 
standing of the resulting halo occupation function, and as 
such it is complementary to the semi-analytical modelling 
described above. 

Clearly, P{N\M) is closely related to how galaxies form 
in dark matter haloes, and so is an important quantity in the 
theory of galaxy formation. However, P{N\M) only contains 
information about the number of galaxies per halo (brighter 
than a certain lower limit). Probably the the most important 
characteristic of a galaxy is its luminosity L. In this paper 
we therefore extend the halo-occupation modeling technique 
by labelling the galaxies with luminosities. We use the 'con- 
ditional luminosity function', $(L|M)dL, which gives the 
average number of galaxies with luminosities in the range 
L ± dL/2 as a function of halo mass M, to address the 
clustering and abundance of galaxies both as function of lu- 
minosity. It is easy to see that the conditional luminosity 
function provides a direct link between the galaxy luminos- 
ity function (hereafter LF) ${L)dL, which gives the num- 
ber of galaxies per comoving volume with luminosities in 
the range L ± dL/2, and the halo mass function n{M)dM, 
which gives the number of dark matter haloes per comoving 
volume with masses in the range M ± dAf/2, according to 

/•oo 

'l'(L) = / $(LiM)n(M)dM. (1) 
Jo 

We use this equation to place constraints on <l>(L|Af) by 
using the halo mass function given by currently favored cos- 
mological models and the observed galaxy LF from the 2dF- 
GRS. It is easy to see that there is an infinite amount of 
different <E>(L|M) that result in exactly the same ^{L) for 
given n{M). We thus need additional constraints, for which 
we use the luminosity dependence of the galaxy-galaxy two- 
point correlation function fgg(r') obtained from the 2dFGRS 
by Norberg et al. (2001a). In addition, we use ^{L\M)dL to 
compute the average total luminosity of galaxies in a halo 
of mass M, 

/>oo 

{L){M)= I $(L|M)LdL (2) 
Jo 

and compare the resulting average mass-to-light ratios as 
function of halo mass, {M/L){M), with observational esti- 
mates on the scales of galaxy clusters and individual galaxies 
(from the TuUy-Fisher relation). This allows us to break the 
model degeneracy and to place interesting constraints both 
on galaxy formation and on cosmology. 

This paper is organized as follows. In Section |^ we first 
give an overview of the observational constraints used in 
this paper. Section |^ presents our model for the conditional 
luminosity function, and the method used to compute lumi- 
nosity functions, mass-to-light ratios and correlation lengths 
as function of luminosity. Results for the ACDM concor- 
dance model are presented in Section ^. In Section ^ we 
discuss possible systematic errors in current observational 
data that may affect our results. Section ^ presents models 
for different cosmologies and for different assumptions about 
the nature of the dark matter. A detailed comparison with 
results obtained from semi-analytical models for the forma- 
tion of galaxies is presented in Section In Section we use 
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the conditional LF to present several statistics regarding the 
distribution of light in the local Universe, and we summarize 
our results in Section ^. 

Throughout this paper we define M to be the halo mass 
inside the radius 7?vir inside of which the average density is 
180 times the cosmic mean density. All luminosities are in 
the photometric fej-band and we convert all units to ft = 1, 
where h = ffo/(100 kms"^). 



2 OBSERVATIONAL CONSTRAINTS 

2.1 The galaxy luminosity function 

Galaxy LFs are typically obtained from large redshift sur- 
veys, the size and accuracy of which have improved steadily 
over the years, (e.g., Efstathiou, Ellis & Peterson 1988; Love- 
day et al. 1992; Marzke, Huchra & Geller 1994; Zucca et 
al. 1997; Folkes et al. 1999). The most recent determinations 
of the LF of "field" galaxies are those from the Las Cam- 
panas Redshift Survey (LCRS; Lin et al. 1996), from the 
SDSS commissioning data (Blanton et al. 2001), and from 
the 2dFGRS (Norberg et al. 2001b). In this paper we use the 
LF of the 2dFGRS, which is currently the largest published 
redshift survey, including more than 110500 galaxies with 
17.0 < bj < 19.2 and z < 0.25. After applying fc-corrections, 
correcting the magnitudes for galaxy evolution and zero- 
point offsets, and correcting the measured number densities 
for incompleteness, Norberg et al. derive a bj LF spanning 
the magnitude range —22.5 ^ h'h- — 51og/i ^ —14.0 which 
is well fit by a Schechter (1976) function 

$(L)dL = |1 (A) " eM-L/L") AL, (3) 

with L* = 9.64 x 10^ /i^"^ L0 (where we have used MQ^b- = 
5.3), a = -1.21 and $* = 1.61 x lO'^/i^Mpc-^ The deriva- 
tion of a luminosity function is cosmology dependent, and 
the parameters listed above correspond to Qq ~ 0.3 and 
= 0.7 0. 

2.2 The galaxy correlation function 

Numerous observational studies in the past have investi- 
gated differential galaxy clustering as function of luminos- 
ity. Although several of these have claimed a dependence 
of clustering strength on luminosity (e.g., Borner, Mo & 
Zhou 1989; Park et al. 1994; Benoist et al. 1996; Willmer 
et al. 1998; Guzzo et al. 2000), others found results that in- 
dicated no significant luminosity dependence (e.g., Loveday 
et al. 1995; Szapudi et al. 2000; Hawkins et al. 2001). The 
main reason for this disagreement is the relatively small size 
of the redshift surveys used. Recently, Norberg et al. (2001a) 
investigated the luminosity dependence of galaxy cluster- 
ing using the 2dFGRS, which is the largest extant redshift 
survey. After correcting their luminosities for band-shifting 
(fe-correction) and galaxy evolution they computed the two- 
point correlation function ^gg{r) for various bins in absolute 
magnitude. Each of these correlation functions is well fit by 

t We also use LFs derived for other cosmological models, which 
were kindly made available to us by Shaun Cole and Peder Nor- 
berg (private communications). 



a simple power-law ^gg{r) = (ro/r)''. Although the shape 
of the correlation function varies little with luminosity (i.e., 
7 ~ 1.7±0.1 for all bins), Norberg et al. find a clear increase 
of the correlation length ro with increasing luminosity. The 
correlation length ro{L) changes slowly for L ^ L* , but in- 
creases by about a factor two from L = L* to L ~ 51/*. 

Note that the conversion of redshift into comoving dis- 
tance is cosmology dependent. The ro{L) and (,gg{r) that we 
use here, and which we adopt from Norberg et al. (2001a), 
are computed assuming flo — 0.3 and JIa ~ 0.7. Note that, 
although the galaxy luminosities have been corrected for 
evolution, the correlation function itself may also evolve with 
redshift, an effect that has not been corrected for. As we 
show below, for the redshift range covered by the 2dFGRS, 
this effect is not negligible. 

2.3 Tully-Fisher relation 

The Tully-Fisher (TF) relation for spiral galaxies relates a 
luminosity measure to a dynamical mass measure and there- 
fore provides an observational constraint on the mass-to- 
light ratio for this class of galaxies. The TF relation can be 
written as 

log(L) = OTP -I- &TP log Vmax (4) 

where Vmax corresponds to the maximum of the observed 
rotation curve and utf and &tp are two constants. This can 
be recast in an expression relating the total mass-to-light 
ratio to the total mass: 

^°^{t) " ~ ^) + (1-797 6tp - OTp) 

-6tp log (-y^) - ^ log^^o (5) 

The third term on the right hand side of equation ^ relates 
the observed rotation measure to the virial velocity Vvir of 
the dark matter halo. Here Vvii is the circular velocity at the 
virial radius within which the average overdensity is 180. 

We use the data obtained and collected by TuUy & 
Pierce (2000; hereafter TPOO). These authors list photom- 
etry and HI linewidths for 24 nearby galaxies for which 
distances are available from Cepheids (the 'local calibra- 
tor' sample), as well as similar data for a large number 
of galaxies in various clusters (the 'cluster' sample). Dis- 
tances to these clusters are obtained by assuming that the 
TF relation is a universal relation. For a sample of in total 
115 galaxies TPOO list the absolute B-band magnitudes and 
the HI linewidths. Note that these magnitudes have been 
corrected for inclination- and linewidth-dependent internal 
dust extinction. Since the luminosities used in the LF func- 
tion and correlation lengths described above have not been 
corrected for internal extinction, we should also use a TF 
relation that has not been corrected for the internal dust 
obscuration. To do this, we use the total apparent _B-band 
magnitudes and distance moduli listed in TPOO and com- 
pute the absolute -B-band magnitudes corrected for galactic 
foreground extinction (based on the 100 /^m cirrus maps of 
Schlegel, Finkbeiner & Davis 1998). These are converted to 
bj-hand magnitudes using bj — B — 0.28 {B — V) (Blair & 
Gilmore 1982) and adopting B — V — 0.7, which corresponds 
roughly to the average color of disk galaxies (de Jong 1996). 
Finally we introduce the scaling with i/o adopting h — 0.7, 
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Figure 1. The left panel plots the fej-band TuUy-Fisher (TF) relation for the 'local calibrator sample' (solid circles) and the 'cluster 
sample' (open circles) uncorrected for intrinsic extinction (data taken from TuUy & Pierce 2000). The solid line indicates the best-fit 
linear regression logL;,^- = 3.79 + 2.671ogVmax- The panel on the right plots the resulting mass-to-light ratios as function of halo mass, 
where we have adopted Vmax/Vvir = 1-6. In addition, the parallelogram indicates the average mass-to-light ratios of disk galaxies obtained 
from the TF relation for six different values of Vmax/Vvir as indicated. See text for a more detailed discussion. 



which corresponds to the value of the Hubble constant used 
throughout this paper. 

The resulting &j-band TF relation is shown in the left 
panel of Figure |l|, where we have made the assumption that 
the HI linewidth is equal to twice the maximum circular ve- 
locity, Knax- Using a double regression with errors in both 
magnitude and ymax, we obtain the best-fit TF relation 
Mt- — 51og^ — —4.17 — 6.671og'Knax, for which otf = 3.79 
and &TF ~ 2.67. These are the values we substitute in equa- 
tion (^) to constrain the mass-to-light ratios of our models. 

For a dark matter halo with a NFW density distri- 
bution (Navarro, Frenk & White 1997) the circular veloc- 
ity increases from the center, reaches a maximum Vinax at 
r ~ 2.163r3 (where Va is the scale radius in the NFW profile) 
and then decreases again to reach Vvii at the virial radius 
Rvii = CTs- One finds that 



14, 



Kir 



: 0.465 



c(l + c) 



(l + c)ln(l-|-c)-c 



(6) 



With our definition of the virial radius (see Section 0) , typi- 
cal galaxy sized haloes have 10 ^ c ^ 50 (Jing 2000; Bullock 

et al. 2001b), corresponding to 1.2 <, Vmax/Kir ^ l-^ The 
maximum rotation velocity of a disk in a halo is typically 
larger than the Knax of the dark matter halo, because of the 
contribution of the baryons to the total circular velocity and 
the contraction of the dark matter halo during the forma- 
tion of the disk (see e.g., Barnes & White 1984; Blumenthal 
et al. 1986; Flores et al. 1993; Mo, Mao & White 1998; Mo 
& Mao 2000; van den Bosch 2002). 

i Note that the c used here is not the same as that in Navarro, 
Frenk & White (1997), Jing (2000) or Bullock et al. (2001b), 
because of our different definition of i?vir- When computing the 
range of expected c we use the results of Bullock et al. (2001b) 
properly scaled to our definition of the virial radius. 



In the right panel of Figure ^, we plot the mass-to-light 
ratios for the individual galaxies of the TPOO data set (open 
and solid circles) computed using Vmax/Vvir = 1-6. In addi- 
tion, we plot the relations between M/L and M obtained 
from the best-fit linear TF relation for 6 different values of 
Knax/Vvir as indicated. We only plot these relations over 
the range of masses that correspond to the magnitude range 
over which the TF relation of TPOO has been derived. In 
what follows we use these TF constraints to judge the cred- 
ibility of our models, keeping in mind that for typical CDM 
cosmogonies one expects 1.4 ^ Vmax/Vvir ^ 2.0 (taking the 
'boosting' effect of the baryons into account). 



3 POPULATING DARK HALOES WITH 
GALAXIES 

3.1 The dark matter mass function 



The mass function of dark matter haloes at z 
written in the form 

diner 



n{M)AM=^vf{u) 



dlnAf 



AM, 



can be 



(7) 



where p is the mean matter density of the Universe at 2 = 
0,1^ = (5c/cr(M), (5c is the critical overdensity required for 
collapse at z = 0, and f(y) is a function of v to be specified 
below. The quantity a{M) in the above equation is the linear 
rms mass fiuctuation on mass scale M, which is given by the 
linear power spectrum of density perturbations P{k) as 



^ ^''^ - 2^ 



P{k) Wiiik) dfc. 



(8) 



where Wjif (fc) is the Fourier transform of the smoothing filter 
on mass scale M . 
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According to the Press-Schechter formalism (Press & 
Schechter 1974), the function /(i^) has the universal form, 



= 2 ( - 



1/2 



(9) 



independent of cosmology, redshift, power spectrum and its 
normalization (e.g.. Bond et al. 1991; Lacey & Cole 1993). 
However, various studies have shown that the halo mass 
function with this /(z^) is inconsistent with numerical sim- 
ulations (e.g., Efstathiou, Ellis & Peterson 1988; Jain & 
Bertschinger 1994; Tormen 1998; Gross et al. 1998; Gov- 
ernato et al. 1999; Jenkins et al. 2001). Using ellipsoidal 
rather than spherical collapse conditions, Sheth, Mo & Tor- 
men (2001) derived an improved form for f{iy) given by 



(10) 



with I/' = y/Eu, a = 0.707, q = 0.3 and A ^ 0.322. 
The resulting mass function has been shown to be in excel- 
lent agreement with numerical simulations, as long as halo 
masses are defined as the masses inside a sphere with an 
average overdensity of about 180 (Sheth & Tormen 1999; 
Jenkins et al. 2001). Therefore, in what follows we consis- 
tently use that definition of halo mass, and we adopt a mass 
function with f(u) given by equation (|l^). In addition we 
use the CDM power spectrum of Bardeen et al. (1986), and 
we adopt a spatial top-hat filter for which 



WM{k;R) 



[sm{kR) - kRcos{kR)] 



(11) 



where the mass M and filter radius R are related according 
to M = 47rpi?^/3. 



3.2 The conditional luminosity function 

In order to compute a LF from the halo mass function, 
we need to specify the conditional luminosity function 
$(L|M)dI/ (see equation Q), which gives the expected num- 
ber of galaxies with luminosities in the range L ± dL/2 (in 
some chosen photometric band) in a halo of mass M. Note 
that $(I/|A/) is a statistical function, and should not be in- 
terpreted as the LF of galaxies residing in any individual 
dark matter halo. 

For massive haloes, such as clusters of galaxies, that 
contain many galaxies the shape of $(L|M) should be the 
same as that of the cluster LF, which can be well described 
by a Schechter function. For all mass haloes, the average of 
$(I/jA/) over the halo mass function should give the field 
galaxy LF. Therefore, we assume that $(LjA/)dL can be 
described by a Schechter function: 



$(LjAf)dL 



£1 

I* 



exp(-L/Z*)dL 



(12) 



Here L* = L*{M), a = a{M) and $* = $*(Af); i.e., the 
three parameters that describe the conditional LF depend 
on M. In what follows we do not explicitly write this mass 
dependence, but consider it understood that quantities with 
a tilde are functions of M. 

With <I>(L|M) defined by equation (|l|), the total aver- 
age luminosity in a halo of mass Al is 
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Figure 2. A comparison of the galaxy LF with the halo mass 
function. Open circles with errorbars correspond to the 2dFGRS 
LF in the fej-band. Solid (dashed) lines correspond to the LF that 
one would obtain from the Sheth & Tormen (Press-Schechter) 
halo mass function under the assumption that each halo yields 
exactly one galaxy with M/L = lOO/i (M/L)0. Note that un- 
der such naive assumptions one expects too many both faint and 
bright galaxies, suggesting that in reality the M/L decreases (in- 
creases) with mass at the low (high) mass end. 



{L){M)= $(L|Af)LdL = l'*Z*r(Q + 2) 



(13) 



with r{x) the Gamma function. The average number of 
galaxies brighter than L* in a halo of mass M is 

$(LlM)dL = 'l*r(5-f-l,l). (14) 

with T{a,x) the incomplete Gamma function. 

For each halo, we define a 'central' galaxy whose lumi- 
nosity we denote by Lc. We assume the central galaxy to be 
the brightest one in a halo, consistent with the fact that in 
most (if not all) haloes the brightest members reside near 
the center. The luminosity of this central galaxy is defined 



Lc(M) 



^{L\M)LdL = -I* L* r{a + 2, Li/L*), (15) 



with Li defined so that a halo of mass M has on average 
one galaxy with L> L\, i.e.. 



$(L|A/)dL : 



(16) 



Note that with this definition, the luminosities of the cen- 
tral galaxies in haloes of the same mass are all the same. 
We have also experimented with drawing Lc at random from 
the distribution function <^{L\M) at L > Li, but found that 
this did not significantly influence our results. We therefore 
adopt the definition of equation (^^ throughout. This is fur- 
ther motivated by numerical simulations and semi-analytical 
models of galaxy formation, which suggest that the lumi- 
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nosities of central galaxies are tightly correlated with the 
masses of their host haloes (e.g. Katz, Weinberg & Hernquist 
1996; Fardal et al. 2001; Pearce et al. 2000; Kay et al. 2001; 
Kauffmann et al. 1993; Somerville & Primack 1999; Cole et 
al. 2000). 

In order to fully specify the conditional LF, we need 
to specify the mass dependence of $*, L* and a. We are 
guided by a direct comparison of the halo mass function 
n{M) and the galaxy LF, ^{L). Under the assumption 
that each dark matter halo harbors exactly one galaxy, 
and that each system has exactly the same mass-to-light 
ratio M/L, the galaxy LF follows directly from the halo 
mass function. Figure Q compares the LF thus obtained for 
M/L = lOO/i (M/L)q with that of the 2dFGRS. The ac- 
tual LF is steeper (shallower) than the one obtained directly 
from the halo mass function at high (low) luminosities. This 
immediately suggests that in reality, rather than being con- 
stant, {M/L) decreases with increasing mass, reaches a min- 
imum at around L* , and then increases again. We therefore 
consider the following simple parameterization. 



q(M) = ai5 log(Afi5). 



(20) 



M\ ,,,, 1 /A/\ \ f M\-l^ ( M 



(17) 



which has four free parameters: a characteristic mass M\ , for 
which the mass-to-light ratio is equal to {M/L)o, and two 
slopes, 13 and 71, which specify the behavior of {M/L) at the 
low and high mass ends, respectively. For M <C Mi one has 
that {M/L) oc M~'^, while for M > Mi the mass-to-hght 
ratio increases with mass according to {M/L) oc M"*"^ (both 
P and 71 are assumed to, but not restricted to, be positive). 

In addition we introduce a fifth free parameter, Mmin, 
which sets the mass scale below which no galaxies form; 
i.e., we set {L){M) = for M < Mmin. This additional 
free parameter has a physical motivation (see below), and 
is included to take account of the fact that the 2dFGRS LF 
has a faint magnitude cut-off at Alt^ — 51ogft — —14.0. 

It is obvious that for galaxy sized haloes L* must 
decrease with decreasing mass, since galaxies with L = 
IO^^/i-^Lq should not be found in a halo ofM = lO^h'^MQ. 
Furthermore, since the LF of the total galaxy population 
has an exponential tail at L > L*, it is unlikely that L* 
keeps increasing with M, because otherwise we would ob- 
serve galaxies brighter than, say, lOOL*. Therefore we adopt 
the following parameterization for L* (M): 



M 



L*{M) 
with 



r(a + 2) 
r(fi + i,i)- 



(18) 



(19) 



This parameterization of L* (M) has two additional free pa- 
rameters: a characteristic mass M2 and a power-law slope 
72. The factor f{&) is to ensure that J\f*{M) = 1 for 
M ^ min[A/i,M2], which is motivated by the fact that 
for small mass haloes (i.e., smaller than a typical galaxy 
group) one expects on average only one dominant galaxy per 
halo. For M 2> max[Ai"i,Af2] this parameterization yields 
Af* oc Afi'2-^i. 

Finally, for 3(Af) we adopt a simple linear scaling with 
log(Af), which we parameterize as follows: 



Here M15 is the halo mass in units of 10^^ h~^MQ, 0:15 = 
a{Mi5 = 1), and rj describes the change of the faint-end 
slope a with halo mass. Note that once a{M) and L*{M) 
are given, the normalization of the conditional LF, $*(Af) 
can be specified through the mass-to-light ratio using equa- 
tions (0) and (^. 

Our model for $(L|A/) specified above thus has 9 free 
parameters: {M/L)o, Afmin, Afi, A/2, /3, 7i) 72, V and 0:15. 
Clearly, it is infeasible to explore the entire 9-dimensional 
parameter space. However, we can use observational con- 
straints and theoretical considerations to significantly re- 
duce the freedom of the model. First of all, we use con- 
straints obtained from studies of the LFs of clusters of galax- 
ies. Beijersbergen et al. (2002) derived the B-band LF of 
the Coma cluster, which is well fit by a Schechter func- 
tion with Z/coma = 6.73 X 10^h~^ Lq and cvcoma ~ —1.32. 
Therefore, we set Q15 = —1.32, which reduces the num- 
ber of free parameters by one. Secondly, I/coma is very sim- 
ilar to L* of the entire galaxy population (see Section 
In addition, the luminosity functions of galaxies in groups 
and poor clusters also have characteristic luminosities sim- 
ilar to L* (e.g. Muriel, Valotto & Lambas 1998). This sug- 
gests that L* does not increase much with mass over the 
ranee 10" - 10^^ h-^MQ, and we therefore tune the value 
of A/2 so that Z*(A/l) = L* with Ml = lO"'^ /i^^Mq. 
Based on a variety of observations, Fukugita, Hogan & Pee- 
bles (1998) find that clusters of galaxies, which they de- 
fine as systems with M > 10^'^h~^ M©, have a mass-to- 
light ratio {M/L)b = 450 ± lOO/i {M/L)q. Therefore, we 
set {M/L){M) = 500/1 (M/L)q for M > lO^/i"^ Mq, and 
tune 71 so that {M / L){M) reaches this upper limit at A/ = 
10"/i^^ Mq. Finally we adopt Afmin = 1.0 x lO^h'^ Mq, 
which corresponds roughly to the mass scale below which 
no galaxies can form in a re-ionized medium, i.e., the virial 
temperature is Tvir ^ 10*K for such systems at the present 
time. 

All in all we have thus reduced the number of fr ee p a- 
rameters to five: (M/L)o, Mi, {3, 72, and r). In Section [l.2| we 
address the impact that this reduction of the model freedom 
has on our results. 



3.3 The galaxy correlation function 

If the initial density distribution of matter is described by 
a Gaussian random field, the clustering properties of dark 
matter haloes are fully determined by the linear power spec- 
trum P(k) and cosmology (Kaiser 1984; Bardeen et al. 1986; 
Cole & Kaiser 1989; Mo & White 1996; Mo, Jing & White 
1997; Sheth, Mo & Tormen 2001). Once a model is adopted 
for the galaxy occupation numbers of dark matter haloes 
(i.e., once P{N\M) is specified), the clustering properties of 
galaxies, such as the two-point correlation function ^gg(r), 
can be computed. With our model for the conditional lumi- 
nosity function <^{L\M) described above, we basically have 
completely specified the differential halo occupation num- 
bers as function of luminosity and halo mass. Therefore, 
once we adopt an additional model for the spatial distribu- 
tion of galaxies in individual haloes, we can compute the 
correlation function of galaxies as function of luminosity. 
Using the definition that 5gg(r-o) = 1 we compute ro{L), 
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which we co mpa re to the observed relation of the 2dFGRS 
(see Section [2.2[ ) . 

We proceed as foUows. We decompose the two-point 
galaxy-galaxy correlation function into two parts: 



(21) 



Here represents the correlation due to pairs of galaxies 
within the same halo (the "1-halo" term), while ^gg describes 
the correlation due to galaxies that occupy different haloes 
(the "2-halo" term). For the 1-halo term we need to specify 
the distribution of galaxies in individual haloes, while the 
2-halo term is given by the correlation of the population of 
dark matter haloes. 



3.3.1 The 2-halo term 

Consider all galaxies with luminosities in the range [Li,L2], 
and define ^gg(r) as the two-point correlation function of this 
subset of galaxies. By definition, the total number of pairs of 
galaxies per comoving volume in this luminosity range that 
have separations in the range r ± dr/2 is 



[1 + Cgg(0]47rr2dr. 



(22) 



where the factor 1/2 corrects for double counting of each 
pair. If ere ng is the mean number density of galaxies with 
L £ [Li, L2] which is given by 



and 



(iV(M)) 



n(M) {N{M)) dM , 



$(L|M) dL 



(23) 



(24) 



gives the mean number of galaxies in the specified luminosity 
range for haloes of mass Af. 

Consider a galaxy with L G [Li,L2] that lives in a 
halo of mass Mi. There are, on average, n(M2) 47rr^dr dark 
matter haloes with mass A^2 that are separated from this 
galaxy by a distance in the range r ± dr/2. In each of these 
haloes there are on average {N{M2)) galaxies that also have 
a luminosity in the range [Li,L2]. Integrating over the halo 
mass function, and taking account of the clustering prop- 
erties of dark matter haloes, we obtain the excess number 
A'^^- of galaxy pairs per galaxy in a halo of mass Mi with 
separations in the range r ± dr/2: 



N+{Mi) = 47rr^dr 

J 

n(M2) {N{M2)) Chh(r; Mi, M2) dMa. 



(25) 



Here 5hh(r;Mi,M2) is the cross correlation between dark 
matter haloes of mass Mi and M2. Note that we implicitly 
assume that all galaxies in a given halo are located at the 
halo center. Since the separation between individual haloes 
is typically much larger than the separation between galaxies 
in the same halo, this simplification does not influence our 
results for large r. 

Multiplying iV+(Mi) with {N{Mi)) and integrating 
over the halo mass function gives the total excess number 
of galaxy pairs per comoving volume with separations in the 
required range. Combining this with equation (M) yields 



1000 



n — I I I 1 1 1 1 1 — 1 — I I I I 1 1 1 1 — I — I I I 1 1 1 1 1 — 3 




0.1 1 

r (h~*Mpc) 



Figure 3. The non-linear dark matter mass correlation func- 
tion $dm(r) of the ACDM "concordance" model (solid line). The 
dashed and dot-dashed lines indicate the "1-halo" and "2-halo" 
terms, respectively. Note that the correlation length ro, defined 
as 5dm (ro) = 1 (dotted line), is determined by the "2-halo" term 
only. 



'^il'^Ar)4nr^dr- 



n{Mi){N{Mi)) N+{Mi) dAfi.(26) 



What remains is to specify the dark halo correlation 
function 5hh(r; Mi, M2). Mo & White (1996) developed a 
model, based on the Press-Schechter formalism, that de- 
scribes the bias of dark matter haloes of mass M with re- 
spect to the dark matter mass distribution. This model was 
extended by Sheth, Mo & Tormen (2001) to account for el- 
lipsoidal collapse. According to this. 



Chh(r; Ml, Af2) = b{Ah) 6(M2) ^Z{r), 



(27) 



where ^dm(^) is the 2-halo term of the dark matter mass 
correlation function (see below) and 



fe(Af) 



(a!/2)= + fo(l-c)(l-c/2) 



(28) 



with a = 0.707, b = 0.5, c = 0.6 and u = 5c/o{M) (Sheth, 
Mo & Tormen 2001). This model has been shown to accu- 
rately match the correlation function of dark matter haloes 
in Af-body simulations (Jing 1998; Sheth & Tormen 1999). 
Substituting equations (^^ and (^ ) in (|2^) we obtain 



^gg('') = 6 idm{r), 
where 



(29) 



b=—l n{M){N{M))b{M)dM . (30) 
% Jo 

Note that the dark matter mass correlation function in the 
above expression is the 2-halo term (i.e., the 1-halo term 
is irrelevant for the halo-halo correlation). Using the same 
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concept as in equation (|2^) we can write that Cdm('') = 
Cdm(»') — 5dm('')- Here ^dm('') is the non- linear two-point cor- 
relation function of the dark matter distribution, and Cc[m(?') 
is the 1-halo term. The Fourier transform of the one-halo 
term is given by 



(fc) = / n{M) [6{M;k)YdM, 



where 

S{M;k)= r^e-*-d^r 
Jn P 



(31) 



(32) 



is the Fourier transform of the halo density profile p(r) (e.g. 
Neyman & Scott 1952; McCleUand & Silk 1977; Ma & Fry 
2000). We assume that p(r) has the NFW form 



p(r) 



Sp 



(r/rs)(l + r/rs 



(33) 



where Ts is a characteristic radius, p is the average density of 
the Universe, and 5 is a dimensionless amplitude which can 
be expressed in terms of the halo concentration parameter 
c — Rvir/rs as 



180 

~3~ ln(l + c)-c/(l-Fc)' 



(34) 



Numerical simulations show that c is correlated with halo 
mass, and we use the relation given by Bullock et al. (2001b), 
converted to the c appropriate for our definition of Rvii. 

The 1-halo and 2-halo terms of the dark matter two- 
point correlation function thus obtained are shown in Fig- 
ure I for a ACDM model with Oq = 0.3, = 0.7, h = 0.7 
and (jg — 0.9. Ifere we have used the fitting formula for 
Cdm('") given by Peacock & Dodds (1996). Note that for 
r ^ 3 h~^Mpc, ^dm{r) is dominated by the 2-halo term, 
while for r ^ 1 /i~^Mpc the 1-halo term dominates. There- 
fore as long as 5 is sufficiently large, ro is determined by the 
2-halo term only. In this case, the value of ro depends only 
on the mean occupation number {N{M)) and is independent 
of the details of how galaxies are distributed in individual 
haloes. In practice, we always find that at r = ro the 2- 
halo term dominates, and we therefore use equation (|2^) to 
compute the galaxy-galaxy correlation function for galaxies 
with L G [Li,L2] from which we obtain ro{L) by solving 
il^iro) = 1. 

As pointed out in Section bl the correlation functions 
obtained by Norberg et al. (2001a) have not been corrected 
for possible redshift evolution. This implies that the ro{L) 
measurements do not correspond to 2: = 0. In fact, since 
more luminous galaxies can be detected out to higher red- 
shift, the correlation lengths of brighter galaxies correspond 
to galajcy populations with a higher median redshift. Note 
that the luminosities have been corrected to 2; = 0, and 
the LF used therefore corresponds to galaxies at the present 
time. Thus, in order to compare model and observations in 
a consistent way, we calculate the model LF at z = and 
the correlation function at some characteristic redshift of 
the sample in consideration. As an approximation, we esti- 
mate the mean redshift z for each volume-limited sample in 
Norberg et al. (2001a), and calculate an effective bias factor 
bes{z) at this redshift from the bias factor b for z — given 
by equation (p9|). Assuming that the number of galaxies is 
conserved in the redshift interval covered by the 2dFGRS 



{Az ~ 0.3 for the brightest sample, and smaller for fainter 
samples), one can show that these two bias factors are re- 
lated by 



cS(Z 



1 + [Z3(0)/D(z)](6- 



(35) 



(Mo & White 1996, 2002; Fry 1996) where D{z) is the linear 
growth rate at redshift z, related to the linear growth factor 
g{z) by D{z) = g{z) + z). We use this effective bias factor 
and the mass correlation function at the mean redshift 1 
to calculate the correlation length ro for each luminosity 
sample. For all the magnitude samples used here, the effect 
is modest, amounting to ^ 10% change in ro. 



3.3.2 The 1-halo term 

For small separations, where the number of pairs are mostly 
due to galaxies in individual haloes, Cgg(r) depends not only 
on the occupation number of galaxies, but also on how galax- 
ies are distributed in individual haloes. 

To model the galaxy distribution in individual haloes, 
we assume that each halo contains one central galaxy (with 
luminosity Lc , see equation ||l^ ) which is located at the halo 
center. The other galaxies, which we call satellite galaxies, 
are distributed with a number density profile that is identical 
to the density profile of the dark matter (equation fj^]). 
This is clearly a simple assumption, because processes such 
as dynamical friction may lead to some segregation between 
galaxies and dark matter particles. 

If the luminosity of the central galaxy Lc G [Li,L2], 
the total number of pairs between the central galaxy and 
satellite galaxies in a halo of mass M is equal to the total 
number of satellite galaxies in the given luminosity range, 
Ncs = Ns{M); otherwise (i.e. when Lc is not in [Li,L2]) this 
number is zero. The total number of distinct pairs of satel- 
lite galaxies is Nbb = ^Ns{Nb — 1). The distribution of pair 
separations between central and satellite galaxies, which we 
denote by /c(r)dr, has the NFW form. The satellite-satellite 
pairs follow a different distribution of separations, denoted 
by /s(r) dr, which can be determined from the NFW density 
profile. The total number of galaxy pairs with separations 
in the range r ± dr/2 can thus be written as 



iVpai.(M)/(r) 



where 



iVes(M)/e(r)+iVcs(Af)/c(r)^ 



iVpair(M) = -{N{M){N{M) - 1)} 



(36) 



(37) 



is the total number of galaxy pairs regardless of separation, 
and /(r) is determined through equation (^6[). The 1-halo 
term of the galaxy correlation function can then be esti- 
mated from 



^4^(r)4.r^dr : 



dr / n{M){Np,ir{M))f{r)dM, 



(38) 



where {A'^pair(Af)) is the mean number of pairs in haloes 
of mass M (see Berlind & Weinberg 2002 for a similar ap- 
proach). 

Note that (A^pair(Af)) depends not only on the mean 
occupation {N{M)), but also on the second moment. Thus, 
in general, the 1-halo term depends also on the deviation 
of the occupation number from the mean. This deviation 
depends on the details of galaxy formation (e.g. Kauffmann 
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et al. 1999; Benson et al. 2000). As a simple model, we make 
the assumption that N{M) has the probability of A'^ + 1 — 
{N{M)) to take the value N and the probability of {N{M))- 
N to take the value iV + 1, if iV < {N{M)) < iV + 1. In this 
case, the mean number of pairs i^ 



(iVpair) = iV (iV) - - iV (iV + 1) 



(39) 



This particular model for the distribution of halo occupa- 
tion numbers is supported by the semi-analytical models of 
Benson et al. (2000), who found that the halo occupation 
probability distribution is narrower than a Poisson distribu- 
tion with the same mean. In addition, Berlind & Weinberg 
(2002) have shown that distributions narrower than Poisson 
are more successful in yielding power-law correlation func- 
tions. 

We use the 1-halo term only to investigate the shape of 
^gg(r), not to compute ro. Since the 1-halo term depends on, 
and is rather sensitive to, the assumptions about the distri- 
bution of galaxies inside dark matter haloes, we primarily 
focus on r o(L ) and only briefly discuss the shape of ^gg(r) 
in Section 4.3. 



4 RESULTS 

Having specified the model and the observational constraints 
we now proceed as follows. After choosing a cosmologi- 
cal model, which sets the mass function n{M)dM and the 
two-point correlation function ^hh(r) of dark matter haloes, 
we specify the conditional LF $(L|M) and compute the 
LP $(1/) and correlation lengths ro{L) of the model, us- 
ing equations (^ and (p^, respectively. Note that in all 
cases discussed here the correlation lengths, defined through 
^gg(ro) = 1, are completely determined by the 2-halo term 
^gg (r) [i.e., fo{L) is independent of our model for the distri- 
bution of galaxies in individual dark matter haloes]. Using 
Powell's multi-dimensional direction set method (e.g., Press 
et al. 1992), we find the values of {M/L)o, Mi, /3, 72 and rj 
that minimize 



Here 



X{^) 



and 



Nil. 

E 

1=1 



$(LO-'l>(L>) 



X^iro) = ^ 



A<I>(L, 



foiLt) - r(){Li) 



Aro(L, 



(40) 



(41) 



(42) 



with Nij, — 35 and Nr = 8 the number of data points for the 
LF and the correlation lengths, respectively. 

Note that with this definition of x^j we are implicitly 
assigning the same weights to the measurements of ^{L) and 
ro{L). This would be correct if there were no systematics in 



3 Note: Berlind & Weinberg (2002), using the same model, 
adopted (Afp^ir) = — !)■ This is only approximately 

correct and can lead to small errors for individual haloes with 
(N) ^ 5. For the computation of ^^g{r), however, this approxi- 
mation is sufficiently accurate. 



the measurements and if the error properties are the same 
in both cases. It is difficult to judge whether this is true in 
the present data, and so our definition of is somewhat 
arbitrary. For instance, since N,j, > Nr and since the errors 
on $(L) are typically much smaller than the errors on ro{L), 
the x^ minimization routine gives much more weight to fit- 
ting the LF than to fitting the correlation lengths. We will 
therefore also consider different relative weights for the LF 
and correlation length measurements, and only use to as- 
sess the relative goodness-of-fit of different models. Tables 1 
and 2 list the parameters for each model that we discuss in 
the text, together with the corresponding values for X^{^) 
and x^(''o)- 

Before proceeding, it is important to address the as- 
sumptions underlying the method. In addition to the trivial 
implicit assumption that all galaxies live inside dark matter 
haloes, we have assumed that (i) $(L|Af) can be represented 
by a Schechter function, and (ii) ^{L\M) is independent of 
environment; i.e., it is assumed that the galaxy properties 
inside a halo of mass M do not statistically depend on the 
halo's larger scale environment (this assumption only influ- 
ences the estimates of ro{L), not the fitting of the LF). Al- 
though assumption (i) is supported by the fact that the LPs 
of groups and clusters of galaxies are reasonably well fit by 
a Schechter function (e.g., Muriel, Valotto & Lambas 1998; 
Beijersbergen et al. 2002; Trentham & Hodgkin 2002), there 
is no strong motivation for adopting a Schechter form for 
haloes of much lower mass. On the other hand, the Schechter 
function has a reasonable amount of freedom (it is specified 
by three free parameters) , and we have tested in great detail 
that our results are robust against different parameteriza- 
tions of these free parameters. Nevertheless, the validity of 
assumption (i) merits more detailed tests, which we intend 
to address in a future paper. Assumption (ii) is supported by 
the work of Lemson & Kauffmann (1999) who used numeri- 
cal simulations to show that the properties of haloes of given 
mass, such as angular momentum and concentration, do not 
significantly depend on environment, and by the extended 
Press-Schechter formalism (e.g. Lacey & Cole 1993) which 
predicts that the merging history of a halo of a given mass is 
independent of its large-scale environment. This assumption 
can ultimately be tested once better statistical samples of 
galaxies become available. 

4.1 The concordance model 

We first consider the standard ACDM model with Qq = 0.3, 
Oa = 0.7, h = 0.7, T = Qoh = 0.21, and as = 0.9. These 
cosmological parameters are consistent with a wide range 
of observations, and so this cosmology has become known 
as the "concordance" model (e.g., Perlmutter et al. 1999; 
Riess et al. 1998, 2001; Melchiorri et al. 2000; de Bernardis 
et al. 2002). 

The solid curves in Figure ^ show the results for the 
best fit model (Ml) for this concordance cosmology. Here- 
after we refer to this model as our 'fiducial' model, and com- 
pare other models with respect to it. The bottom left panel 
compares the model LF with that of the 2dFGRS (open 
circles with errorbars). The agreement is excellent, which 
is also apparent from the fact that x^{^) divided by the 
number of degrees of freedom is very close to unity. The 
agreement with the observed ro{L) (shown in the bottom 
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Table 1. Model parameters for the ACDM "concordance" model. 





logMi 




"15 


logTV^i 


iOgiW2 




P 


71 


72 


V 




X \^o) 


t 

sgg 






1,3} 






\°) 


(1) 


\°) 


1,9) 


(lOj 




(L2) 


(13j 


(A A\ 


Ml 


13.5 


500 


-1.32 


11.27 


11.73 


134 


0.77 


0.32 


0.65 


—0.50 


41.6 


5.4 




Mia 


13.5 


500 


-1.32 


11.34 


11.73 


133 


0.75 


0.33 


0.65 


—0.50 


41.4 


5.8 


+ 


Mlb 


13.5 


500 


— 1.32 


10.62 


11.77 


145 


0.96 


0.25 


0.66 


—0.53 


90.0 


2.2 


_1_ 

\ 


Mlc 


13.5 


500 


-1.32 


11.88 


11.25 


79 


0.46 


0.52 


0.63 


-0.65 


54.5 


3.0 


+ 


Mid 


13.5 


500 


-1.32 


10.42 


11.74 


102 


0.60 


0.28 


0.69 


-0.36 


40.9 


1.9 


+ 


MIA 


12.5 


500 


-1.32 


8.64 


12.60 


163 


2.81 


0.15 


0.83 


-0.12 


46.1 


26.2 


+ 


MIB 


14.5 


500 


-1.32 


12.82 


12.19 


111 


0.56 


0.80 


0.87 


-0.31 


196.5 


43.9 




MIC 


13.5 


250 


-1.32 


10.37 


12.60 


325 


1.21 


0.05 


0.73 


-0.26 


54.3 


96.1 




MID 


13.5 


750 


-1.32 


11.81 


11.49 


111 


0.57 


0.52 


0.63 


-0.64 


44.2 


3.1 


+ 


MIE 


13.5 


500 


-1.12 


11.06 


11.93 


136 


1.07 


0.29 


0.72 


-0.31 


36.0 


4.7 




MIF 


13.5 


500 


-1.52 


11.46 


11.53 


133 


0.64 


0.34 


0.60 


-0.70 


49.2 


6.6 





Column (1) lists the ID by which we refer to each model in the text. Columns (2) to (4) list model parameters that are kept fixed; here 
Mi is defined such that L*{Ml) = L* , {M/L)^i is the mass-to-light ratio of haloes with M > 10^'^ h~^M.Q, and ais is the faint-end 
slope of the conditional LF for haloes with M = 10^^h~^ Mq. Columns (5) to (11) list the best-fit model parameters, of which M2 and 
71 are set by Ml and {M/L)^i, respectively. Columns (12) and (13) list the values of x^{^) ^i^d X^(^o) of the best-fit model, respectively. 
Finally, column (14) indicates whether the shape of the galaxy correlation functions are in good (-I-), reasonable (-1-/—) or poor (— ) 
agreement with the data. Masses and mass-to-light ratios are in Mq and h (M/L)0, respectively. 



panel on the right) is also acceptable, although the model 
predicts slightly too large correlation lengths for galaxies 
with Mtj — 51og/i J; —20. Note that the functional form of 
ro{L) is remarkably well reproduced by the model. 

The best-fit value of rj is —0.50, which indicates that 
a decreases with increasing mass. Larger values of a im- 
ply a LF that is more peaked at around L* , i.e., with rel- 
atively fewer galaxies with L < L* . Note that a faint-end 
slope that decreases with increasing mass is at least qual- 
itatively consistent with observational results (e.g., Tren- 
tham & Hodgkin 2002). The top left panel of Figure § plots 
the characteristic luminosity, L* , as function of mass. A 
best-fit value for (3 of 0.77 impUes that L* oc M^ '^'^ for 
M < Ml ~ 1.9 X 10"/i"^ M©. The slope of L*(M) for 
haloes with M ^ Mi is mainly controlled by 72, and is 
strongly constrained by the form of ro{L); Larger (smaller) 
values of 72 cause ro{L) to increases less (more) rapid with 
increasing luminosity. Note that for our best-fit value of 
72 = 0.65 the value of L* only increases by a factor ~ 5 
from PI = lO^^ft-i Mq to M = lO^^^/j-i ^g. 

The middle panel on the left plots N{> L*/10), the 
average number of galaxies brighter than one-tenth of the 
characteristic luminosity L* of the 2dFGRS. For M ^ 
10^^ M0 this number increases roughly linearly with 
mass, while below ~ lO^^/i"^ M© one basically no longer 
finds galaxies as bright as L*/10. The top and middle pan- 
els on the right plot the mass-to- light ratios Ad/L and AI/Lc, 
respectively. Here L is the average total luminosity of galax- 
ies in a halo of mass M (see equation |^), and Lc cor- 
responds to the average luminosity of the 'central galaxy' 
(equation [^). The mass-to-light M/L reaches a minimum 
value of ~ 125/1 (M/L)0 at M ~ lO^'^/i"^ Mq. For haloes 
with masses smaller or larger than ~ 10^^'^ M©, the 
galaxy formation efficiency is reduced, consistent with the 
fact that feedback can reduce the star formation efficiency 
in small haloes, while cooling becomes less effective in more 
massive haloes (e.g.. White & Rees 1978; Dekel & Silk 
1986; van den Bosch 2002; Benson et al. 2000). Note that, 



by construction, the mass-to-light ratio stays consta nt a t 
M/L = 500/i(M/L)q for M > 10^*/i"^Mq (see Section |3l|) . 

For low mass haloes, the ratio M/Lc has properties sim- 
ilar to M/L. This owes to the fact that, by construction, each 
of these haloes hosts only one dominant ('central') galaxy. 
For massive systems, the faint-end slope of the conditional 
LF is steeper and M* > 1, causing M/Lc > M/L. The 
parallelogram indicates the TF constraints for six different 
values of Vmax/Vvir from 1.0 (upper line) to 2.0 in steps of 
0.2 (see Figure |l|). In order for model Ml to be consistent 
with the zero-point of the TF relation, the ratio Vmax/K'ir 
must be smaller than 1.4, which is only marginally consistent 
with CDM predictions, especially when the boosting effect 
of the baryonic matter is taken into account (see discussion 
in Section 2.3). In addition, the fact that M/Lc{M) is quite 
strongly 'curved' over the mass range of the TF constraints 
implies that Knax/K'ir needs to vary significantly, and in a 
fine-tuned fashion, with mass in order to yield a linear TF 
relation^]. Thus although model Ml yields a reasonable fit 
to the LF and ro{L), it predicts values for (M/Lc) that are 
slightly too high. 

As pointed out before, by virtue of the definition of 
used in the minimization routine, much more weight 
is given to fitting the LF than to fitting the correlation 
lengths. Consequently, the fits are expected to be controlled 
mainly by the LF constraints. To demonstrate this, the dot- 
ted curves in Figure ^, corresponding to model Mia, show 
the results obtained when fitting only the LF, i.e., by setting 
~ x^(*&)- The fact that these curves are all very similar 
to our fiducial model (solid curves) confirms the above as- 
sertion. In addition, the dashed curves in Figure ^ corre- 
sponding to model Mlb, show the results obtained by set- 
ting x^ ~ X^(^)/50 + X^i^o), i.e. by reducing the weights 



'I Note that the discrepancy for haloes with M ^ lO^^h-i Mq 
may not be a too serious problem, because their central galaxies 
may be predominantly ellipticals for which the TF relation does 
not apply. 
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Figure 4. Results for the fiducial model Ml (solid lines). The two bottom panels compare the model predictions with the observed 
LF (left) and correlation lengths ro(L) (right). Open circles with errorbars correspond to the 2dFGRS data from Norberg et 

al. (2001a, b). The four top panels show, clockwise from the upper left panel, the characteristic luminosity, L^., the ratio between halo 
mass and total luminosity M/L, the ratio between halo mass and luminosity of the 'central' galaxy M/Lc, and the number of galaxies 
with L > L,/10, all as function of halo mass. The dotted and dashed curves show the results obtained by minimizing = x^('J') (Model 
Mia) and = X^(*)/50 + x^(^o) (Model Mlb), respectively. 



of the measurements by a factor 50. Not surprisingly, 

this model fits the correlation lengths better but fits the 
LF slightly worse, at least according to the value of x^(*&)- 
Whether or not the fit is acceptable depends on possible 
systematic errors in the observed $(-£/) which may be larger 
than the statistical errors used here. We return to the issue 
of systematic errors in Section 5.1. 



4.2 Less restrictive models 

The results presented above suggest that, for the cosmo- 
logical concordance model, the conditional luminosity dis- 
tribution obtained from the observed LF and the observed 
clustering strength as function of luminosity over-predicts 
the mass-to-light ratios implied by the TuUy-Fisher rela- 
tion. However, this assertion is only correct for our restricted 
set of parameterized models for $(L|M). The more general 
model for our conditional LF has nine, rather than five, free 
parameters. Although we have argued that four of these pa- 
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Figure 5. Results of models in which one model parameter is changed relative to the fiducial choice. Left panels show M/Lc as 
function of halo mass (parallelogram indicates the TF constraints), while the right panels show the correlation length as function of 
absolute magnitude. Open circles with errorbars correspond to the observational data from the 2dFGRS. Solid curves in each panel 
correspond to the fiducial model Ml. Top panels show results for models with = 10^^'^ h~^'MQ (Model MIA, dotted curves) and 
Ml = Mq (Model MIB, dashed curves). Middle panels show results for models with (M/L)ci = 250h (M/L)© (Model MIC, 

dotted curves) and (M/L)ci = 750h (M/L)© (Model MID, dashed curves). Bottom panels show results for ais = -1.12 (Model MIE, 
dotted curves) and ais = —1.52 (Model MIF, dashed curves). See text for a detailed discussion. 



rameters can be fixed based on physical and observational 
constraints, we now address the impact these 'restrictions' 
have on our results. 

In our fiducial model we fix the mass M2 by set- 
ting Ml = 10^^-^h~^ Mq, where Ml is defined so that 
L*{Ml) ~ L* . This choice is motivated by the fact that L* 
does not seem to vary much over the range 10^^ /i~^M0 ^ 
M <, 10^^ /i"^M0. To illustrate the impact that this as- 



sumption has on our results, the upper panels of Figure g 
plot the results for models MIA and MIB where we set 
Ml = W^^-^ /i"^M0 and Ml = lO"'^ h'^MQ, respectively. 
Note that increasing (decreasing) the value of Ml mainly 
causes bright galaxies to reside in less (more) massive haloes. 
As one can see, the results depend quite sensitively on the 
value of Ml- For Ml > lO^^'^ H'^Mq the predicted ro{L) 
only shows a very mild luminosity dependence, in clear dis- 
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agreement with the data. In addition, this model yields val- 
ues of X^(^) that are too high and predicts galaxy-galaxy 
correlation functions that deviate strongly from the power- 
law form for Mi, — 51og/i ^ —21.0 (see Section 4.3 below). 
Model MIA, for which Ml = lO^^ '^ fe-^M©, also fails to 
match the observed roi^L) for bright galaxies. Both models 
predict mass-to- light ratios that are clearly inconsistent with 
the TF constraints. We therefore conclude that models with 
10" ft"^M0 <, Ml^ 10" '' /i"^M0 are the most successful, 
and in what follows we keep Ml fixed at lO" '' h'^MQ. 

In the fiducial model, we tune 71 so that M/ L — 
(M/L),i = 500/1 (M/L)q at M = lO^^h-^ Mq, and we keep 
M/L fixed at the value of {M/L)ci for more massive systems. 
The middle two panels of Figure ^ show the effects of chang- 
ing (M/L)ci to 250/i(M/L)q (model MIC) and 75O/i(M/L)0 
(model MID). These values are still consistent with the ob- 
servational constraints of Fukugita et al. (1998) at the So- 
level. Increasing (decreasing) (M/L)ci causes a decrease (in- 
crease) of the number of galaxies in massive haloes. Since the 
clustering strength of haloes increases with mass, an increase 
(decrease) of (Af/L)ci therefore results in lower (higher) val- 
ues of ro{L). For (M/L) a < 500/i (M/L)© the models pre- 
dicts too large correlation lengths as well as mass-to-light 
ratios that are inconsistent with the TF constraints. Model 
MID, however, predicts a ro{L) and "I>(L) that are both in 
agreement with the data. Over the mass range probed by 
the TF constraints, the mass-to-light ratios are almost iden- 
tical to those of the fiducial model; they are slightly too high 
and reveal an uncomfortably large variation over this mass 
range. Since (M/L)ci = 750/i (M/L)© is only marginally 
consistent with the observational constraints of Fukugita et 
al. (1998), we restrict ourselves to (M/L) a = 500/i (M/L)© 
in what follows. 

The lower panels of Figure ^ compare the fiducial model 
(for which ais = -1.32) to models MIE and MIF for which 
Q15 = —1.12 and ais = —1.52, respectively. The choice of 
Q15 has only a small impact on the results. A change in 
ai5 is largely compensated by changes in Mi, /3 and 77, but 
does not significantly improve or worsen the agreement with 
the data. H owev er, 015 does impact on the shape of Cgg(^) 
(see Section 4.3 below) with ais = —1.32 clearly giving the 
best results. We therefore restrict ourselves to models with 
Q15 = —1.32 in what follows. 

Finally we investigated the effect of changing Mmin, i.e., 
the mass scale below which "1>(L|M) = 0. We find that as 
long as Mmin <, 10" ft-^M© our results do not significantly 
depend on Mmin, and our fiducial value of lO^/i~^M0, which 
is motivated by re-ionization considerations, therefore does 
not infiuence any of the results presented here. 

In addition to varying some of the parameters kept fixed 
in our fiducial models, we also experimented with completely 
different parameterized forms of L* , a, and Neither of 
these alternative models, however, was able to yield results 
in better agreement with the data. Therefore, we conclude 
that the ACDM concordance model can be made consis- 
tent with both the observed LF and the correlation lengths 
as function of luminosity. The largest problem for the cos- 
mological concordance model is to fit the zero-point of the 
TF relation. In all the models presented here, the mass- 
to-light ratios are only consistent with the TF constraints 
if Vmax/Vvir Si 1.4, which is hard to reconcil e w ith typi- 
cal CDM predictions (see discussion in Section 2.3). Proba- 



bly equally problematic, is the fact that the models predict 
mass-to-light ratios that vary strongly over the mass range 
probed by the TF constraints. 

4.3 The shape of the correlation function 

It is well known that the observed correlation function of 
galaxies is well fit by a single power law ^gg(r) = (r/ro) ' 
over a wide range of radii (e.g. Davis & Peebles 1983; Baugh 
1996; Jing, Mo & Borner 1998; Zehavi et al. 2002). Norberg 
et al. (2001a) have recently shown that even when ^gg{r) 
is split in several luminosity bins, each individual corre- 
lation function is still well fit by a single power law with 
7 ~ 1.7. These observations yield extra constraints on the 
conditional luminosity function. Note that the mass correla- 
tion function of the CDM model in consideration is not well 
described by a pure power law (see Fig.^. The power-law 
behavior of ^gg(r) for individual luminosity bins therefore 
results from some non-trivial conspiracy between the halo 
bias on large scales (the 2-halo term) and the distribution 
of galaxies within individual dark matter haloes (the 1-halo 
term) on small scales. This property of the galaxy occupa- 
tion model has been noticed in some early investigations 
(e.g. Jing, Mo & Borner 1998; Jing, Borner & Suto 2002; 
Kauffmann et al. 1999; Benson et al. 2000; Berhnd & Wein- 
berg 2002). 

Figure ^ plots the galaxy-galaxy correlation functions of 
model Ml (solid lines) for five luminosity bins also used by 
Norberg et al. (2001a). In addition, we also plot predictions 
for the magnitude bin — 24 < Mt- — 5\ogh < —23. For all 
sub-samples with Mb — 51ogh, > —22.5, the correlation func- 
tions have roughly a power-law form over the entire range 
from r ~ 0.1h~^ Mpc to r ~ 20h~^ Mpc, and with a slope 
close to that measured from the 2dFGRS (indicated in each 
panel). For galaxies with —23.0 > Mtj — 5\ogh > —24.0, 
however, the correlation function is strongly suppressed at 
r <, 2/i~^Mpc. These bright galaxies correspond to the cen- 
tral galaxies of rich clusters, and so their correlation function 
contains only a negligible contribution from the 1-halo term. 
Unfortunately, the 2dFGRS results of Norberg et al. (2001a) 
only considered galaxies with Mbj — 51og/i > —22.5, so that 
we cannot compare this prediction with data. Figure ^ also 
plots the predictions for models MIE (0:15 = —1.12) and 
MIF (ai5 = —1.52). Both these models show significant 
deviations from a single power-law form (at least for the 
brighter galaxies), demonstrating that the actual shape of 
the correlation function can be used to discriminate between 
different models for the conditional LF. It is remarkable that 
the model with ais = —1.32, which is the value obtained by 
Beijersbergen et al. (2002) for the Coma cluster, yields the 
best result. 

Although the shapes of the two-point correlation func- 
tions of our fiducial model are in good agreement with the 
2DFGRS, this success is model-dependent. In particular, we 
have assumed that the distribution of galaxies in individ- 
ual haloes follows the density distribution of dark matter, 
and we have made simplified assumptions regarding the sec- 
ond moment of the distribution of the halo occupation num- 
bers. Indeed, as shown in Berlind & Weinberg (2002), the 
shape of the correlation function on small scales may de- 
pend sensitively on these assumptions, and a good power 
law is not always obtained. Therefore, in judging the suc- 
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Figure 6. The predicted two-point correlation functions of galaxies in different absolute magnitude bins. Results are shown for the 
fiducial model (Model Ml; solid curves), and for models MIE (ais = —1.12; dotted curves) and MIF (ais = —1.52; dashed curves). 
The short solid line-segments indicate the slopes of the best-fit power-laws obtained by Norberg et al. (2001a) for 2dFGRS galaxies in 
the same magnitude bins. No such line-segment is shown in the lower right panel, since no 2dFGRS data is available for this magnitude 
bin. Note that the correlation functions of the fiducial model resemble power-laws over the entire radial range plotted, with slopes that 
are very similar to those observed. An exception is the magnitude bin in the lower right panel, which reveals a clear deviation from a 
power law at small separations. Models MIE and MIF, on the other hand, show significant deviations from a pure power law behavior 
for galaxies with -21 > M;,^ - 51og/i > -22.5. 



cess of our models we will not give too much weight to the 
actual shape of (,gg{r). For completeness, we do indicate in 
Tables 1 and 2 whether the correlation functions, obtained 
under the assumptions of our model for the 1-halo term, are 
in good reasonable (+/— ) or poor (— ) agreement with 
the data. 



5 SYSTEMATIC ERRORS IN THE 
OBSERVATIONAL CONSTRAINTS 

As we have shown above, none of the ACDM concordance 
models is fully consistent with all observational constraints. 
Here we address whether this discrepancy might result from 
uncertainties and/or systematic errors in the observational 
constraints. 
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Figure 7. Same as Figure |4| except that here we compare the fiducial model Ml with models in which we assume that the true LF 
differs from the 2dFGRS LF. In particular, we show results for models in which we assume that the true number density of galaxies is 
1.5 times as high as for the 2dFGRS LF, independent of luminosity (Model Mlc; dotted curves), and in which we assume that the true 
faint end slope is a = —1.4 in stead of —1.21 (Model Mid; dashed curves). Note that both models results in lower values of M/Lc (in 
better agreement with the TF constraints) and in smaller correlation lengths (in better agreement with the 2dFGRS data). 



5.1 The luminosity function 

The determination of an accurate LF is non-trivial. The ac- 
curacy depends, amongst others, on the size of the redshift 
survey (the larger the survey, the smaller the systematic 
errors due to large scale structure), on surface brightness 
completeness effects, on the accuracy of star-galaxy separa- 
tion, and on the accuracy of the photometric calibration. In 
addition, corrections have to be made for band-shifting (the 
fc-correction) and for evolutionary effects of the galaxy popu- 
lation with redshift, each of which has its own uncertainties. 



Our main motivations for using the 2dFGRS LF are that 
(i) it is the largest extant redshift survey, (ii) its quoted er- 
rors have contributions from all the effects mentioned above, 
and, most importantly, (iii) the differential galaxy cluster- 
ing properties, used as additional constraints on the model, 
have been derived from the same 2dFGRS data set. 

Despite these advantages over other LFs available in the 
literature, one cannot rule out that some amount of system- 
atic errors are present. For instance, Blanton et al. (2001) 
claim that the 2dFGRS LF (initially presented by Folkes et 
al. 2001) is inconsistent with the LF of the SDSS commis- 
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sioning data. Norberg et al. (2001b), however, claim that 
after taking account of zero-point photometric offsets, after 
correcting for galaxy evolution, and after making corrections 
for the proper conversion between SDSS and 2dFGRS mag- 
nitudes, both LFs are actually in good agreement. 

Despite this apparent concordance, it is worthwhile to 
investigate how the presence of certain systematic errors im- 
pacts on our results. One of the most likely sources for sys- 
tematic errors is incompleteness. Note that the incomplete- 
ness corrections applied by Norberg et al. were based on a 
comparison with the SDSS early-data release, which itself is 
not expected to be 100 percent complete. If indeed the 2dF- 
GRS suffers from incompleteness effect, the true luminosity 
densities will be higher, resulting in lower mass-to-hght ra- 
tios for individual dark matter haloes. This is exactly the 
effect that might alleviate the problems with the TF con- 
straints outlined above. We therefore address the impact 
incompleteness may have on our results using two different 
approaches. 

First we consider an overall, luminosity independent in- 
completeness, by assuming that <1>* of the 2dFGRS is a factor 
1.5 too small (Model Mlc). Note that this is roughly similar 
to the disagreement between the SDSS and 2dFGRS LFs 
claimed by Blanton et al. (2001). Secondly we consider the 
possibility of a luminosity-dependent incompleteness by as- 
suming that the true value of a is —1.4 rather than —1.21 
(Model Mid). Results are shown in Figure ^ Not surpris- 
ingly, both models result in mass-to-light ratios in better 
agreement with the TF relation. Thus, the ACDM concor- 
dance model can be made consistent with all observational 
constraints considered here if the 2dFGRS is only about 70 
percent complete at all luminosities. Given the discussion in 
Norberg et al. (2001b) such extreme incompleteness effects, 
however, seem unlikely. 

The incompleteness effects discussed here have exactly 
the same effect as the presence of dark CDM haloes (i.e., 

haloes in which no observable galaxies formed)|j||. In this 
case, model Mlc is to be interpreted as a model in which 
one third of all haloes (independent of mass) host no ob- 
servable galaxies (i.e., "1>(L|M) — 0), while in the case of 
model Mid the fraction of dark CDM haloes increases with 
decreasing halo mass. The possibility of dark CDM haloes 
was recently discussed by Verde, Oh & Jimenez (2002), who 
argued that haloes with large angular momenta might result 
in disks of such low surface brightness that star formation is 
never triggered (see also Jimenez et al. 1997). As shown by 
Jimenez, Verde & Oh (2002) a large abundance of dark CDM 
haloes is consistent with the rotation curves of disk galaxies. 
Based on our discussion above, we therefore conclude that 
the disagreement of the ACDM concordance model with the 
TF constraints may be alleviated if a relatively large frac- 
tion of dark matter haloes never managed to form observable 
galaxies. 



5.2 The correlation lengths 

As we have argued above, the differential clustering strength 
of galaxies as function of luminosity provides important ad- 



Note that we have already assumed that haloes with M < 



ditional constraints that allow us to obtain a unique best-fit 
$(I/|M). In particular, the value of 72 is strongly constrained 
by the strong increase in correlation length at the bright end. 

In this paper, we used the results from the 2dFGRS, 
presented by Norberg et al. (2001a), for the reasons already 
mentioned in Section 5.1. A comparison of the ro{L) ob- 



tained from the 2dFGRS with those obtained from a variety 
of other surveys, in particular the Southern Sky Redshift 
Survey 2 (Benoist et al. 1996; Willmer et al. 1998), the 
Stromlo/APM redshift survey (Loveday et al. 1995), and 
the ESO Slice Project (Guzzo et al. 2000), shows an overall 
agreement that ro{L) increases with increasing luminosity 
(see Figure 3 in Norberg et al. 2001a). Quantitatively, how- 
ever, the disagreement is relatively large compared to the 
quoted errorbars. Undoubtedly, to a large extent this owes 
to sampling fluctuations associated with the relatively small 
size of some of these surveys (e.g., Benson et al. 2001). Since 
the 2dFGRS is currently by far the largest redshift survey for 
which differential clustering measures have been obtained, 
we adhere to these data. In the near future, when the full 
2dFGRS and SDSS results are available, the true accuracy 
of the ro{L) measurements used here can be addressed in 
more detail. 



5.3 The Tully-Fisher zero-point 

We have argued that the AI/Lc of our fiducial model are at 
best marginally consistent with the (zero-point of) the TF 
relation. We computed our TF constraints from the data 
of TPOO, for which the absolute magnitudes have been de- 
termined using distances obtained with the Cepheid period- 
luminosity relationship. The zero-point of this TF therefore 
intrinsically contains an effective Hubble parameter. In or- 
der to be able to make a comparison with our models and 
with the data on the LF and the galaxy-galaxy clustering, 
we converted the absolute magnitudes of the TPOO data to 
Mb^ ~ 51og/i adopting Ho = 70 kms"^ Mpc~^ This partic- 
ular choice for the Hubble constant is motivated by (i) the 
value Ho — 72 ± 8kms~^ Mpc~^ (68% confidence) obtained 
by the HST Key Project (Freedman et al. 2001), (n) the 
value Ho = 77± 8kms~^ Mpc~^ (95% confidence) obtained 
by TPOO from the same TF data used here, and (iii) the fact 
that we have adopted h — 0.7 in all our models. 

However, some uncertainty regarding Ho remains, and 
this will directly influence the normalization of the TF zero- 
point. For example, determinations of Ho based on time- 
delay measurements of gravitational lenses usually yield val- 
ues for the Hubble parameter that are signiflcantly smaller 
(e.g. Kochanek 2002; but compare Hjorth et al. 2002). 
Kochanek (2002) obtains Ho = 48+4 kms'^Mpc"^ (95% 
confidence) under the assumption that the lenses have flat 
rotation curves. Furthermore, Saha et al. (2001), using 
Cepheid calibrations for a sample of galaxies with Type la 
supernovae, finds Ho = 59 ± 6 kms^^ Mpc^^. 

It is straightforward to correct our TF constraints for 
an error in Ho- The mass-to- light ratios obtained from equa- 
tion (H) scale as 



'-(f)„--(f)„-'-(^) 



(43) 



Thus, for h < 0.7 the mass-to-light ratios obtained from the 
TPOO TF relation become larger. In other words, if for a 
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Table 2. Model parameters for different cosmologies. 



ID 






^8 


n 




mx 


logA/i 


logA/2 


(A//L)o 


13 


71 


72 


V 


x'(*) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


M2 


0.2 


0.8 


1.06 


1.0 


- 


- 


11.25 


11.18 


72 


0.60 


0.42 


0.64 


-0.73 


60.5 


72.4 


+ 


M3 


0.3 


0.7 


0.86 


1.0 


- 


- 


11.42 


11.69 


132 


0.72 


0.34 


0.65 


-0.52 


41.6 


2.5 


+ 


M4 


0.4 


0.6 


0.74 


1.0 






11.10 


12.16 


226 


1.02 


0.22 


0.66 


-0.40 


45.8 


12.9 


+ 


M5 


0.5 


0.5 


0.66 


1.0 






10.92 


12.58 


369 


1.40 


0.14 


0.69 


-0.31 


45.9 


42.9 


+ 


M6 


1.0 


0.0 


0.46 


1.0 






8.74 


14.34 


5674 


2.89 


-0.14 


0.83 


-0.09 


305.4 


187.0 


+ 


M7 


0.3 


0.7 


0.80 


1.0 






11.66 


11.62 


127 


0.64 


0.38 


0.64 


-0.55 


42.3 


1.1 


+ 


M8 


0.3 


0.7 


0.70 


1.0 






12.40 


11.49 


120 


0.46 


0.57 


0.61 


-0.61 


45.7 


3.5 


+ 


M9 


0.3 


0.7 


0.90 


0.8 






11.32 


11.57 


110 


0.70 


0.36 


0.65 


-0.55 


44.1 


28.2 


+ 


MIO 


0.3 


0.7 


0.90 


1.2 






11.20 


11.88 


160 


0.84 


0.28 


0.66 


-0.46 


41.2 


3.3 


+ 


Wl 


0.3 


0.7 


0.90 


1.0 


0.1 


0.8 


10.93 


11.64 


118 


0.64 


0.30 


0.66 


-0.52 


40.7 


8.4 


+ 


W2 


0.3 


0.7 


0.90 


1.0 


0.2 


0.5 


10.00 


11.37 


70 


0.25 


0.29 


0.68 


-0.55 


48.4 


10.2 


+/- 



Column (1) lists the ID by which we refer to each model in the text. Models whose ID starts with M (W) correspond to CDM (WDM) 
models. Columns (2) to (5) indicate the parameters of the cosmological model, where the values modified with respect to the fiducial 
model (the ACDM concordance model) are in bold face. Columns (6) and (7) indicate the free streaming scale (in Mpc) and the 
dark matter particle mass (in keV), respectively (for WDM models only). Columns (8) to (17) are the same as columns (5) to (14) in 
Table 1. 



given halo mass AI and our assumption that h — 0.7 the 
M/Lc of a model implies a value (Vmax/K'ir)o.7, then for 
h ^ 0.7 the implied ratio is 

iK-Jr (K-)o..no:7) 

where we have used &tf = 2.67 (see Section [2. j ). For exam- 
ple, in our fiducial model (Ml) the values of M/Lc imply 
that Knax/K/ir ^ 1.4, which is hard to reconcile with the 
typical circular velocity curves of CDM haloes. If, however, 
h — 0.5 then we obtain V^ax/Vvir ^ 1-8, in much better 
agreement with CDM predictions. 



6 COSMOLOGICAL CONSTRAINTS 

So far we have only focussed on the ACDM concordance cos- 
mology. In the present section, we turn our attention to a 
number of alternative cosmologies to (i) investigate whether 
the method presented here can be used to constrain cosmo- 
logical parameters and (ii) to see whether we can find a set 
of cosmological parameters that improves the fit to $(!/) 
and ro(L) with respect to the concordance model while si- 
multaneously yielding mass-to-light ratios that are in better 
agreement with the TF constraints. 

In a recent study, Zheng et al. (2002) have shown that 
distinct cosmological models produce distinct halo popula- 
tions. Based on a detailed investigation of halo populations 
in a wide range of cosmologies, Zheng et al. argued that a 
combination of galaxy clustering and statistics that are sen- 
sitive to the virial masses of dark matter haloes should be 
able to impose strong constraints on cosmological models 
without having to rely on a priori assumptions regarding 
galaxy formation. Indeed, recently Guzik & Seljak (2002) 
and Yang et al. (2002) used the galaxy-galaxy weak lensing 
resuhs of the SDSS (McKay et al. 2001) to constrain the 
virial masses of galaxies brighter than L* (see also McKay 
et al. 2002). Linking these results to the abundance of dark 
matter haloes, Seljak (2002) was able to place useful con- 
straints on cosmological parameters. This suggests that our 



combination of constraints from the luminosity-dependent 
clustering strengths and the mass-to-light ratios implied by 
the TF relation, should impose constraints on cosmological 
parameters. As we show below, this is indeed the case. 



6.1 The matter density 

We consider cosmologies in which the matter density is var- 
ied with respect to the concordance value Qq — 0.3 (models 
M2 - M6; see Table 2). All models considered here have 
SIa ~ 1.0 — ^0, to keep agreement with the cosmic mi- 
crowave background constraints (e.g., Melchiorri et al. 2000; 
de Bernardis et al. 2002), and we adopt F = Qoh throughout. 
In addition, when we modify flo we also change the normal- 
ization as, in order to satisfy the weak-lensing constraints 
obtained by Hoekstra, Yee & Gladders (2002): 

as^a^^ = 0.46 . (45) 

This allows us to investigate whether our modeling of the 
conditional LF can break the degeneracy between as and 
fio that results from weak lensing and cluster abundance 
measurements. 

Comparing models with Qo ^ 0.3 with data is some- 
what complicated by the fact that obtaining $(-L) and 
ro(-L) from observations is cosmology dependent. The ob- 
servational constraints used so far were obtained assuming 
Slo = 0.3 and f^A = 0.7, consistent with the cosmological 
parameters considered up to this point. For other values of 
fio and Qa, however, the observational constraints from the 
2dFGRS wiU be diflterent. For the LF, S. Cole and P. Nor- 
berg kindly provided $(1/) obtained from the 2dFGRS for 
the additional cosmologies considered here (private commu- 
nications). For the ro{L), however, only measurements for 
fio = 0.3 and Qa = 0.7 are available to us (from Norberg et 
al. 2001a). We therefore ignore the cosmology dependence of 
the observed correlation lengths in our analyses. Since the 
Slo dependence is small compared to the statistical errors on 
ro(L), this does not significantly affect our results. Finally, 
the TF constraints also depend on Qq (see equation |^). This 
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Figure 8. Same as Figure ^ except that here we compare results for a number of ACDM models with different values for £7o as indicated. 
All models have Ha = 1.0 — Oq, normalization o-sOq'®^ = 0.46 (Hoekstra et al. 2002), and shape parameter F = floh with h = 0.7 (see 
Table 2). Note that the mass-to-light ratios (upper two right panels) are plotted as function of M' = M (flo/O.S)^'^^ , and that the 
middle panel on the right plots log[M'/ic] rather than log[M/Lc]. This takes care of the fact that the TF constraints depend on 
(see equation |fe|) and allows all models to be plotted in the same Figure. Note also how model with Qq > 0.35 yield mass-to-light ratios 
that are too high, while models with Qq < 0.35 result in correlation lengths that are too large. 



effect has been taken into account in the middle right panel 
of Figure ^ by rescaling M/Lc and M accordingly so that 
we can plot the curves for different flo in the same figure. 

Figure ^ presents results for the best fit models with 
f2o = 0.2, 0.3, 0.4, 0.5 and 1.0. Models with Qq < 0.3 pre- 
dict values of ro{L) that are too large, while models with 
Qo ^ 0.3 predict too high values of M/Lc to match the ob- 
served TF constraints. Apparently, the technique used here 
is extremely effective in breaking the degeneracy between Qo 
and erg that is present in current measurements of cosmic 



shear and/or the abundance of rich galaxy clusters. Remark- 
ably, the best results are obtained for ilo close to the con- 
cordance value of fio = 0.3, in good agreement with the con- 
straints obtained from high redshift supernova la searches 
(e.g., Riess et al. 1998; Perlmutter et al. 1999) and from 
observations of the temperature anisotropies in the cosmic 
microwave background (e.g., Sievers et al. 2002). 

For an Einstein-de Sitter cosmology, we find mass-to- 
light ratios that are more than an order of magnitude too 
high to match the TF constraints. This is in good agreement 
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log[M] (h-^Mg) Mb-51ogh 

Figure 9. Same as Figure ^ except that here we compare models in which we vary as (upper panels), the slope of the initial power 
spectrum n (middle panels), and the nature of the dark matter (lower panels) with respect to the fiducial model. All other cosmological 
parameters are kept the same as for the ACDM concordance cosmology. 



with previous studies based on a variety of different methods 
(e.g.. White & Frenk 1991; Cole 1991; Kauffmann et al. 1993; 
Steinmetz & Navarro 1999; Elizondo et al. 1999; van den 
Bosch 2000). 

6.2 The normalization of the power spectrum 

In the previous section we have shown that models with 
ilo = 0.3 yield the best results, once erg is kept normalized 
to the constraints obtained by Hoekstra et al. (2002). For the 
ACDM concordance model this implies erg = 0.86. This nor- 
malization is consistent with constraints from the observed 



abundance of rich clusters of galaxies (e.g. White, Efstathiou 
& Frenk 1993; Eke, Cole & Frenk 1996; Bahcall & Fan 1998; 
Viana & Liddle 1996, 1999; PierpaoU, Scott & White 2001) 
and from weak gravitational lensing studies (van Waerbeke 
et al. 2002; Bacon et al. 2002; Refregier et al. 2002). How- 
ever, a number of recent studies have suggested a somewhat 
lower normalization, of the order of erg = 0.7 - 0.8 (e.g., 
Carlberg et al. 1997; Viana, Nichol & Liddle 2002; Seljak 
2001; Reiprich & Borhringer 2002). We therefore investigate 
the impact of such small changes in erg, while keeping the 
value of 0,0 fixed at 0.3. 

Changes in erg mainly affects the clustering strength, 
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which is proportional to 6^cr|. The bias b{M) increases with 
decreasing a{M) (Mo & White 1996), and since a{M) oc 
erg, lowering the value of as causes haloes of given mass to 
be more strongly biased. However, since on galactic scales 
the increase of b with l/a{M) is less rapid than linear, the 
clustering strength decreases when as is lowered. The upper 
panels of Figure ^ compare our fiducial model Ml (erg — 0.9) 
with models M7 {as = 0.8) and M8 (erg = 0.7). Indeed, lower 
values of erg result in smaller correlation lengths. In fact, for 
as ~ 0.8 the model yields an excellent fit to the observed 
ro{L)- In addition, from the values of X^(^) it clear that 
these modest changes of ag do not significantly influence 
the quality of the fit to the LF. However, the M/L ratios 
on the mass scale probed by the TF relation become slightly 
higher, therewith only worsening the problem of reproducing 
the correct TF zero-point. Thus, lowering erg by about 10 
percent with respect to the concordance value brings the 
luminosity dependence of the galaxy clustering strengths in 
excellent agreement with observations, but an alternative or 
additional modification is required to bring the mass-to- light 
ratios in agreement with the TF zero-point. 



6.3 Tilted CDM models 

So far we have only considered models with a Harrison- 
Zel'dovich power spectrum (n = 1). However, inflationary 
models allow small deviations of n from unity, and we there- 
fore investigate the effect of small modifications of n. 

The middle panels of Figure ^ plot the results for mod- 
els M9 (n = 0.8) and MIO (n = 1.2). In both models we keep 
the normalization of the power spectrum fixed at ag = 0.9 
and we adopt the cosmological parameters of the concor- 
dance model. Lowering n has the effect of suppressing the 
linear power spectrum on small scales, and so one might 
hope to reduce the number density of galactic haloes. We 
find this effect to be quite small, and so the mass-to-light 
ratios required to match the luminosity function are quite 
similar to those obtained for the concordance model with 
n — 1. The biggest effect of changing n is to change the 
correlation length of the mass, for example, the mass corre- 
lation length is about 5.1 h~^Mpc for n = 0.8, and about 
4.7 h~^Mpc for n = 1.2. Since the bias factor for galax- 
ies increases with decreasing n (because of the decrease in 
a{M) on galaxy scales), the predicted correlation lengths for 
galaxies increases with decreasing n, as shown clearly in the 
middle right panel. Thus, an increase of the power index n 
has a similar effect as lowering erg; it brings the correlation 
lengths in better agreement with the data, while slightly 
increasing the mass-to- light ratios on the scale of galaxies. 



6.4 Warm Dark Matter 

In recent years there have been various claims that the CDM 
paradigm is in conffict with observations on small scale. 
First, CDM predicts too many low mass haloes which leads 
to problems with the abundance of satellite galaxies around 
a typical Milky Way sized galaxy (e.g., Moore et al. 1999; 
Klypin et al. 1999; but see Stoehr et al. 2002) and to prob- 
lems with the formation of disk galaxies (White & Navarro 
1993; Navarro & Steinmetz 1999, 2000; Mo & Mao 2000; van 
den Bosch, Burkert & Swaters 2001). Second, CDM models 



predict dark matter haloes with steeply cusped density dis- 
tributions, inconsistent with the observed rotation curves of 
some dwarf and low surface brightness galaxies (e.g., Blais- 
Ouellette, Amram & Carignan 2001; de Blok et al. 2001; but 
see van den Bosch & Swaters 2001). 

In order to (partially) alleviate these problems. Warm 
Dark Matter (WDM) has been suggested as an alternative 
dark matter candidate (e.g., Sommer-Larsen & Dolgov 2001; 
Colin, Avila-Reese & Valenzuela 2000; Bode, Ostriker & 
Turok 2000) . In WDM models the dark matter particles have 
a non-negligible initial velocity, which causes the damping of 
small-scale fluctuations due to free-streaming. Consequently 
the power-spectrum of WDM is modified with respect to 
that of CDM, and can be expressed in terms of the CDM 
power spectrum PcDM(fc) according to 

-PwDM(fc) = TwDM(fc) -PcDM(fc) (46) 

where the WDM transfer function is approximated by 



7wDM(fc) = exp 



fc-R/,WDM (fc-R/,WDM)^ 



(47) 



(Bardeen et al. 1986; Sommer-Larsen & Dolgov 2001). The 
comoving free-streaming scale -R/,wdm is related to a free- 
streaming mass M/_wDM according to 



M 



/,WDM 



Mr n 



"WDM 



-R/.WDM \ 

/i-i Mpc ) 



(48) 



which in turn is related to the mass mx of the WDM par- 
ticles as 



mx = 2Ah-"*a 



5/4ol/2 



M 



/,WDM 



WDM \ IQll/i-lMo 



-1/4 



keV. 



(49) 



The lower two panels of Figure ^ compare the results 
for our fiducial CDM model Ml with two WDM models, Wl 
and W2, that differ only in the mass of the WDM particles 
(see Table 2). Except for the properties of the dark mat- 
ter, these models have the same cosmological parameters as 
the ACDM concordance model. Decreasing mx causes a de- 
crease in the abundance of small mass haloes. Consequently, 
in order to preserve the abundance of (faint) galaxies, each 
(low mass) halo needs to harbor more (or brighter) galaxies. 
This causes a decrease in the mass-to-light ratios. Further- 
more, WDM haloes are expected to be somewhat less con- 
centrated than their CDM counterparts (e.g.. Bode, Ostriker 
& Turok 2000), favoring lower values of Knax/K/ir. Both of 
these two effects bring the model in better agreement with 
the TF constraint. On the other hand, the change from CDM 
to WDM does not significantly change the fit to the lumi- 
nosity function and correlation length. It thus seems that 
WDM models with mx Si 0.8 are consistent with all data 
considered here. However, such low WDM particle masses 
may be hard to reconcile with the observed opacity distri- 
bution of the Lya forest (Narayanan et al. 2000), with the 
phase-space densities of dwarf galaxies (Dalcanton & Hogan 
2001), and with re-ionization constraints (Barkana, Haiman 
& Ostriker 2001). It thus remains to be seen to what ex- 
tent WDM may be considered a successful alternative for 
the CDM paradigm. 
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7 COMPARISON WITH SEMI-ANALYTICAL 
MODELS FOR GALAXY FORMATION 

The problem of simultaneously reproducing the galaxy LF 
and the TP zero-point has been noticed before in studies 
based on semi-analytical models for the formation of galax- 
ies. Initially these studies focussed on Einstein-de Sitter cos- 
mologies, where the discrepancy was found to be very large 
(e.g., Kauffmann et al. 1993; Cole et al. 1994). This is in 
excellent agreement with our results, which indicate that 
cosmologies with fio = 1 imply mass-to-light ratios that 
are more than an order of magnitude too large. Better re- 
sults were obtained for cosmologies with lower values of fio, 
but the overall agreement remained unsatisfactory (Heyl et 
al. 1995). 

Further improvement was achieved by including dust 
obscuration in the modeling (Somerville & Primack 1999; 
Cole et al. 2000; Benson et al. 2000). As long as Kiax/Kir = 
1 was adopted, these models found good agreement with 
the TF relation. However, using the more realistic values in 
accord with the NFW profiles of CDM haloes, the semi- 
analytical models predict too high rotation velocities at 
given luminosity. Similar results were obtained by studies 
that do not attempt to simultaneously reproduce the galaxy 
LF, but normalize the mass-to-light ratios of galaxy-sized 
haloes by adopting baryon fractions motivated by nucle- 
osynthesis constraints (Steinmetz & Navarro 1999; Mo & 
Mao 2000; Navarro & Steinmetz 2000; van den Bosch 2000). 
Once again, our conclusions are in excellent agreement with 
these results. 

In general, the semi-analytical models use the extended 
Press-Schechter formalism (Bond et al. 1991; Bower 1991; 
Kauffmann & White 1991; Lacey & Cole 1993; Somerville 
& Kolatt 1999) to compute Monte-Carlo realizations of the 
halo merger histories. Since this method yields no informa- 
tion about the spatial distribution of dark matter haloes, the 
clustering properties of galaxies can not be addressed. Vari- 
ous studies have therefore applied the semi-analytical mod- 
els directly to dark matter haloes grown in cosmological sim- 
ulations (Kauffmann et al. 1997; 1999; Diaferio et al. 1999; 
Benson et al. 2000, 2001, Springel et al. 2001; Mathis et 
al. 2001) and computed correlation functions of the model 
galaxies. 

Benson et al. (2000) and Kauffmann et al. (1999) pre- 
sented the mass-to-light ratios as function of virial mass 
for their models, which can be directly compared to the 
{M / L){M) obtained from our conditional LF. Benson et 
al. (2000) find a clear minimum of ~ 125/i (M/L)© at 
M « lO^^/i"^ Mq, while for M Z 10"/i"^ Mq their mass- 
to- light ratios converge to a value of ~ 450ft (M/L)©. The 
agreement with our fiducial model (see upper right panel of 
Figure ^ is impressive, suggesting that their semi-analytical 
models result in a conditional LF not too different from the 
one derived here using a completely different approach. This 
is further supported by the fact that their models yield a 
luminosity dependence of the correlation length that is in 
reasonable agreement with our fiducial model. 

This suggests that the technique introduced here, in 
addition to constraining cosmological parameters, can be 
used to place constraints on galaxy formation. For example, 
Kauffmann et al. (1999) showed that different treatments 
of star formation and feedback have a strong infiuence both 



on the slope and the amplitude of the two-point correlation 
function of galaxies. Benson et al. (2000), on the other hand, 
concluded that the predictions for the correlation function 
are robust to changes in the semi-analytical model parame- 
ters as long as the models match the bright end of the lumi- 
nosity function. Our results, however, clearly indicate that 
different models, which fit the LF equally well, can have 
quite different clustering properties. For example, models 
MIA to MIF all match the LF roughly equally well (with 
the possible exception of MIB) but yield correlation lengths 
that differ quite strongly. This suggests that, as concluded 
by Kauffmann et al. (1999), the physics of galaxy formation 
(which set the conditional LF) has a strong effect on the 
clustering properties of galaxies. 

The above discussion clearly indicates the potential 
strength of the technique introduced here. Note that our 
conclusions are independent of the uncertainties present in 
current theories of galaxy formation. In particular, no uncer- 
tain assumptions have to be made regarding the efficiencies 
of star formation and feedback, or regarding the effects of 
merging and dust obscuration. This makes the results pre- 
sented here more robust. An additional advantage is the 
computational costs; only a few minutes are required to 
find the best-fit parameters of the conditional LF for a spe- 
cific cosmological model. Thus, our approach can be used to 
probe a large number of cosmological models in a relatively 
short time. On the other hand, the approach explored here 
has its limitations. It only yields a parameterized form of the 
conditional LF, and it is unclear if the best-fit ^{L\M) can 
be understood in terms of the physical processes governing 
galaxy formation. This may be most easily addressed with 
semi-analytical models, and as such the two approaches are 
complementary. 



8 THE DISTRIBUTION OF LIGHT IN THE 
LOCAL UNIVERSE 

The conditional luminosity function $(L|M) that we have 
attempted to constrain in this paper is a powerful tool in 
linking the distribution of galaxies to that of the dark mat- 
ter. As an illustration, we use the <I>(L|M) of our fiducial 
model to present some statistical properties of the distribu- 
tion of (blue) light (of galaxies) in the local Universe (i.e., at 
z ~ 0). Note that these distributions only give an approxi- 
mate description of the true distribution of light. After all, 
they correspond to our fiducial model, which is not perfectly 
consistent with all data. The description given here therefore 
only serves as an illustration of the statistical distributions 
that can be obtained once the conditional LF is known more 
accurately. 

First we use $(LjAf) to define the joint probability dis- 
tribution 

P{L,M)dLdM = — n(M)$(L!M)LdLdM (50) 

PL 

which gives the probability that a 6j-band photon that is 
produced in the local Universe originates from a galaxy with 
luminosity in the range L±dL/2 that resides in a halo with 
virial mass in the range M ± dM/2. Here pi, is the average 
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Figure 10. Statistics regarding the distribution of blue light at 2 = 0. The upper left panel plots the fraction of blue light emitted by 
galaxies inside haloes with virial masses in the range M, M + dM. The corresponding cumulative distribution is shown in the upper right 
panel, indicating that ~ 50 percent of all blue light is produced inside haloes with M < 2 X 10^^h~^ Mq. The contour plot in the lower 
left panel corresponds to the joint probability distribution P{L, M) (equation [pO||). The cross marks the maximum of P{L, M), and the 
nth contour level is at 2~" times the maximum. Note that galaxies with L ~ lQ^h~^ Lq living inside haloes with Af ~ 1 X 10^^h~^ Mq 
are the most efficient producers of blue light in the local Universe. Finally, the bottom right panel plots the conditional probability 
distributions P{M\L) (equation [p^) for four different luminosities as indicated. 



luminosity density, given by 

/■oo poo 

PL = / dLL dM n{M) <i>{L\M) 
Jo Jo 

= / $(L)LdL. (51) 

The contours in the lower left panel of Figure ^ show the 
joint probability distribution P{L, M) for our fiducial model 
Ml. As can be seen, the joint probability distribution peaks 
at around L ~ L* and M ~ 1 x 10^^ /i^^M©, i.e., these sys- 
tems are the most efflcient producers of blue light in the local 



Universe. Integrating P{L, M) over all luminosities yields 
the fraction of all light, fL^M) AM , produced by galaxies 
inside haloes with virial masses in the range M ± AM/ 2: 

fL{M)= [ P{L,M)AL= —n{M)-{L){M) (52) 
Jo PL 

The upper left panel of Figure |l^ plots as function of halo 
mass again for the conditional LF of our fiducial model. Note 
that haloes with M ~ 10^^ h~^MQ are the most produc- 
tive in 'producing' light. This is not surprising, given that 
M / L(M) reaches a minimum at around this mass scale. This 
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is also consistent with generic models of galaxy formation, 
which predict that haloes more massive than lO^^ft-^Mg are 
inefficient in cooling their baryonic gas, while less massive 
haloes are inefficient in producing stars because of heating 
from supernova feedback (e.g., Larson f974; White & Rees 
1978; Dekel & Silk 1986). 

The upper right panel of Figure |l^ plots the cumula- 
tive distribution F{< M) = /l(M) dM, which gives the 
fraction of all light produced in haloes less massive than 
M. In our fiducial model, about 50 percent of all light is 
produced inside haloes with M ^ 2 x 10^^ h~^MQ. Since 
this is roughly the mass scale where one changes from iso- 
lated haloes around galaxies to haloes that harbor groups 
and clusters of galaxies, this indicates that isolated galax- 
ies produce roughly the same amount of light as galaxies 
residing in clusters and groups. Note also that the distri- 
bution of /_L(Af) is quite narrow. Galaxies in haloes with 
M ^ 2 X 10^* /i'^Mq and M ^ 5 x 10^° h'^MQ each are 
responsible for only about 5 percent of the total blue light 
in the local Universe. 

Finally, we use the conditional LF to compute the con- 
ditional probability distribution P{M\L)dM that a galaxy 
of luminosity L lives inside a halo with virial mass in the 
range M ± AM/2. Using Bayes' theorem, 

P{M\L) AM = ^^^^ n{M) AM . (53) 

The lower right panel of Figure |l^ plots this probability 
distribution for four different luminosities: L — L*/100, 
L = Z/*/10, L = L* , L = 10 L*. Whereas lOL* galaxies are 
typically found in haloes with lO^^/i"^ ^ M <, lO^'^^"^ Mq, 
galaxies with L = L'/lOO ~ 10^h~^ Lq typically reside in 
haloes of M ~ 5 x 10^" h~^MQ. Only a very small fraction 
of these galaxies {^5 percent) are found as a member of a 
larger group or cluster of galaxies. It should be pointed out, 
however, that our model predictions for P{M\L) are quite 
uncertain for galaxies with L ^ L*/10. This owes to the 
fact that constraints on ro{L) only exist for galaxies with 
L > L* /5, even though data on the luminosity function ex- 
tends down to much fainter luminosities. If, for example, 
for L < L*/5 the correlation lengths increase with decreas- 
ing luminosity (i.e., opposite to the trend for the brighter 
galaxies), we would obtain conditional luminosity functions 
with many more faint galaxies in massive haloes than in the 
current prediction. Clearly, observational data on the cor- 
relation functions of faint galaxies are needed in order to 
better constrain their P{A4\L). 



9 CONCLUSIONS 

We have presented a novel technique to link the distribution 
of galaxies to that of dark matter haloes. This technique 
extends upon earlier modeling of halo occupation number 
distributions by labelling the galaxies with luminosities. In 
particular, we have introduced the conditional luminosity 
function <&(L|M)dL, which gives the average number of 
galaxies with luminosities in the range L ± AL/2 that re- 
side in a halo of mass M. Starting from the number density 
of CDM haloes, n{M)AM, we sought the $(L|M) that re- 
produces the 2dFGRS fej-band luminosity function $(L)dL. 
Since there are many different <l>(LjM) that are consistent 



with both n(M) and $(!/), we used the luminosity depen- 
dence of the galaxy correlation lengths obtained from the 
2dFGRS as additional constraints. This allows the recovery 
of a unique best-fit <E>(L|A/). To assess the viability of the 
conditional LF thus found, we used $(L|M) to compute the 
average mass-to-light ratios of dark matter haloes as func- 
tion of mass, which we compared to constraints obtained 
from the (zero-point of) the TF relation. 

For the ACDM concordance cosmology {Qq ~ 0.3, 
Qa ~ 0.7, h — 0.7, as = 0.9), we are able to accurately 
reproduce the observed LF. In addition, adopting a sim- 
ple model for the distribution of galaxies inside dark mat- 
ter haloes, we obtain galaxy-galaxy correlation functions as 
function of luminosity that are in agreement with the data 
of Norberg et al. (2001a). However, the values of (AI/L) we 
obtained are slightly too high to match that implied by the 
TF relation. 

We have shown that the concordance model can be 
made consistent with the data if either the 2dFGRS is only 
~ 70 percent complete (which is very unlikely), or if ~ 30 
percent of all dark matter haloes have failed to form observ- 
able galaxies (they only contain 'dark' galaxies). Although 
feedback and an ionizing background can explain the pres- 
ence of such dark galaxies in low mass haloes, there is cur- 
rently no clear mechanism that can prevent star formation 
in more massive haloes, rendering this option unlikely as 
well (but see Verde et al. 2002). 

The mass-to-light ratios of the concordance model can 
also be brought into agreement with the TF constraints if 
Ho is lowered to ~ 55 kms~^ Mpc~^. Although this value 
is only marginally consistent with the recently advocated 
value of i?o = 72 ± 8 kms'^Mpc"^ (68% confidence) of 
the HST Key Program (Freedman et al. 2001), such low 
value is in agreement with recent time-delay measurements 
in gravitational lens systems (Kochanek 2002; but see also 
Hjorth et al. 2002). Alternatively, (M/L) may be lowered 
by adopting warm, rather than cold dark matter. In WDM 
cosmologies the abundance of low mass haloes is suppressed, 
resulting in lower mass-to-light ratios; WDM haloes are also 
expected to be less concentrated than CDM haloes and so 
the mass-to-light ratios implied by the observed TF relation 
is also higher. Both of these effects lead to better agreement 
between the model predictions and the TF constraints. In 
order to have a significant effect, we find that the WDM 
particles must have a mass of the order of mx ^ 0.8keV, 
which is already the lower limit set by a number of other 
observational constraints (Narayanan et al. 2000; Dalcanton 
& Hogan 2001; Barkana et al. 2001). 

In addition to these modifications of the concordance 
cosmology, we have also investigated cosmologies with dif- 
ferent matter densities. We have shown that the technique 
presented here allows one to break the degeneracy between 
Qo and erg inherent in the cosmological constraints obtained 
from cosmic shear and/or cluster abundance measurements. 
The best results are obtained for Qo — 0.3, in good agree- 
ment with the high redshift supernova results (e.g., Riess 
et al. 1998; Perlmutter et al. 1999): for Z 0.4 the mass- 
to-light ratios of dark matter haloes are too high to match 
the TF zero-point, while cosmologies with flo ^ 0.2 yield 
correlation lengths that are too large. 

In summary, the modeling of the conditional LF pre- 
sented here strongly favors a cosmology with Qq ~ 0.3. How- 
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ever, in order to bring the models into agreement with the 
observed zero-point of the TF relation, either h is lowered by 
15 to 20 percent with respect to the concordance values of 
0.7, dark matter is warm rather than cold, or a large fraction 
of haloes host no luminous galaxies. 

Once a satisfactory model for the conditional LF is 
found, it can be used to extract useful statistical informa- 
tion about the distribution of light in the Universe. As an 
example, we have computed the fractional contribution to 
the total luminosity density by galaxies with luminosity L 
that live in haloes of mass M. In particular, we have shown 
that 50 percent of all the blue light is produced in haloes 
with M ^ 2 X 10^^ /i"^ Mq. In addition, the conditional LF 
can be used to compute the probability P{M\L)dM that 
a galaxy with luminosity L lives in a halo with mass in 
the range M ± dM/2. These probability distributions may 
prove useful for a statistical interpretation of the mass-to- 
light ratios inferred from, for example, galaxy kinematics or 
gravitational lensing. 

The technique presented here can easily be extended to 
compute higher-order correlation functions, as well as pair- 
wise peculiar velocity dispersions, all as function of galaxy 
luminosity. In the near future several large redshift sur- 
veys of galaxies will become available, such as the SDSS 
and the completed 2dFGRS at z = and DEEP and VIR- 
MOS a,t z — 1. These will allow far more stringent con- 
straints to be placed on cosmological parameters, especially 
when combined with the complementary constraints from 
weak-lensing and cluster abundances. In addition, once lu- 
minosity functions and correlation lengths are available in 
different photometric passbands and for different morpho- 
logical classes, the technique presented here can easily be 
extended to simultaneously model the distribution of galax- 
ies as function of color and/or morphological type. This will 
not only help to constrain cosmological parameters, but also 
to provide important insights into the complicated processes 
associated with galaxy formation. 
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